^DA0251  90 


RADC-TR-76-112 
Final  Technical  Report 
April  1976 


PREDETECTION  COMPENSATED  IMAGING 
Hughes  Research  Laboratories 


Sponsored  by 

Defense  Advanced  Research  Projects  Agency 
ARPA  Order  No.  2646 


Approved  for  public  release; 
distribution  unlimited. 


./  n 


ROME  AIR  DEVELOPMENT  CENTER 
AIR  FORCE  SYSTEMS  COMMAND 
8RIFFISS  AIR  FORCE  BASE,  NEW  YORK  13441 


SECTIONS  V and  VI  (pages  89  - 135)  have  been  deleted  from  this  document. 

This  report  has  been  reviewed  by  ’_.ie  RADC  Information  Office  (01)  and 
is  releasable  to  the  National  Technical  Information  Service  (NTIS) . At  NTIS 
it  will  be  releasable  to  the  general  public  including  foreign  nations. 

This  report  has  been  reviewed  and  is  approved  for  publication. 


APPROVED: 


DONALD  W.  HANSON 
Project  Engineer 


The  views  and  conclusions  contained  in  this  document  are  those 
of  the  authors  and  should  not  be  interpreted  as  necessarily 
representing  the  official  policies,  either  expressed  or  implied, 
of  the  Defense  Advanced  Research  Projects  Agency  or  the  U.  S. 
Government. 


Do  not  return  *^his  copy.  Retain  or  destroy. 


PREDETECTION  COMPENSATED  IMAGING 


C.  J.  Swigert 
W.  P.  BroWi 
T.  R.  O'Meara 
D.  Y.  Tseng 
S . Wandsura 


Contractor:  Hughes  Research  Laboratories 

Contract  Number:  F30602-75-C-0026 

Rtfecti-ve  Date  of  Contract;  16  September  19/T 

Contract  Expiration  Date:  15  August  1975 

Amount  of  Concract;  5235,000.00 

Program  Code  Number:  5E20 

Period  of  work  covered:  .Sep  74  - Aug  75 

Principal  Investigator:  C.  J.  Swigert 

Piione:  213  456-6411 


Project  Engineer: 

Phone : 


Donald  W.  Hanson 
315  330-3145 


Approved  for  public  release; 
distribution  unlimited. 


This  research  was  supported  by  the  Defense  Advanced 
Research  Projects  Agency  of  the  Department  of 
Defense  and  was  monitored  by  Donald  W.  Hanson  (OCSE) , 
Griffiss  AFB  NY  13441. 


-UNCLASSIFIED 


SECURITY  CL*>^S$‘ Fixation  of  THI'=  n»ta  Emoted) 


4 


f 

t.  REPd^rsUMDER 


^ REPORT  DOCUMENTATION  PAGE 


RAI 


:/tR-76-1]2  / 


4 — Jwgc  f«n»  »u»lfTn«r^ 


L 


U GOVT  ACCeSSJCN  ^0, 


READ  INSTRUCTIONS 
BEFCKE  COMPI.ETING  FORM 


3 5r.O'  ry’”'.  CATAI.OG  NUMBER 


PREDETECTION  COMPENSATED  IMAGING 
- ^ 


1 


• " 

C.  J./Swigert 
W.  P./ Brown  •. 

T.  RJo'Meara 
wti? 


D.  yJ  Tseng 
S.^Wandzura 


— ^ ' 


^ ! W>  Vf/TTttTRyKT  K 0 A OCR  £SS~ 


Fical  technical  Rep«rt» 
16  Sep  74  <— -15  Aug  75  ^ 


/ 

/ 


■« — reRROfisiNj  ORO  wgPcf«IX3>fflTgETr' 

N/A 


8 CONI  RACT  OR  GRANT  NUMSERfsJ 


1 F3MO2-75-C-O026 3 f ^ 

1 ir  > 


wqpi^UNIT  numbers 

Hughes  Research  Laborato^^s  j "*■ 

3011  Malibu  Canyon  RoacT  /'"  Fy  ~ '/  / ^ — ■\'4>6 

Malibu  CA  90265 ^ y ^ j j ''j  i 12 


II.  CONTROuUNG  OFFiCE  NAME  A>7g/rB'^j>689  ' 

Defense  Advanced  Research  Prejects  Agency 
1400  Wilson  Blvd 
Arlington  VA  22209 


l«  MONITORING  AGENCY  NAME  A ACDPESSil/  tfi»/e.-.N,r  .'-o-r  Con.->//-i«  O/i.ee. 

Rome  Air  Development  Center  (OCSE) 

Griffiss  AFB  NY  13441 


2646031 


?PXJ0AI£- 


ApraS3Bg7f  J 


■'Vi'  } / 


180 

15  SSCuftiTv  Class  (of  thie  report) 


UNCLASSIFIED 


IBs.  CSCLASSl'TCATION  DOWNGRADING 
SCHEDULE 


16  OISTRibu'ION  statement  ('n/ (Ala  ReptH 


Approved  for  public  release;  distribution  unlimited 


17  O’STRIBUTION  STATCIMEjajT  (of  tho  in  Block  50,  !i  diffe’ent  from  Report) 

Same 


18.  supplementary  notes 

RADC  Project  Engineer:  Donald  W.  Har.con  (OCGE) 

SECTIONS  V end  VI  (oagvs  89  - 135)  have  been  deleted  from  this  document. 


’5»  KKY  WOftOS  (Continue  on  rev'«r*»  Je  it  '*ec«s' ■ary  and  fdentity  hy'  btock  number) 

COAT,  Predetectiou  Compensation,  Image  Formation,  Turbulence  Compensation, 
Adaptive.  Optical  Systems,  Deforiiiable  Mirror,  Signal  to  Noise  Ratio,  Predictive 
Track,  Visual  Assessment 


20  iBST1?ACT  fCot  *ir>oe  on  revarso  aide  if  neces€mr“  artf*  by  block  nurr.berl 

P'This  report  describes  the  te^hni.jues  developed  and  evaluated  for  predetection 
compensated  imaging  systems  ta  measure  ant  correct  wavefront  distortions  at  a 
telescope  aperture.  One  approach  to  compensate  these  wavefront  distortions  is 
to  introduce  the  appropriate  .'real-time)  wavefront  correction  ahead  of  the 
image  formation  process.  

(Cont'd) 


DD 


1‘!73  EOITION  of  I NCV  65  !S  OR50LETE 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (WhenDtla  Bnlrrrd) 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAGEfHTien  D*(a  Enitrtd) 


We  call  systems  that  achieve  such  corrections  Imaging  COAT  (coherent 
Optical  Adaptive  Techniques)  systems,  since  they  are  spatitfally 
coherent  systems. 

Two  general  classes  of  wavefrent  errer  sensors  have  been  identified 
and  explored  in  this  study.  We  have  also  examined  several  wavefront 
perturbation  techniques  (multidither  techniques)  that  allow  a detector 
to  operate  in  a time  shared  mode  to  collect  the  wavefront  error  data 
that  would  normally  require  the  use  of  several  detectors  using  a non- 
perturbation mode  of  operation.  This  technique  has  been  developed  and 
refined  into  the  Multiple  Aperture  Interlinked  (MAIL)  predetection 
compensation  system. 

When  advantageous,  we  use  the  MAIL  concept  for  predictive  track  of 
multiple  nhase  layers  to  obtain  in  extended  integration  time;  otherwise, 
we  rely  on  *quick  looks*  of  the  phase  error  while  viewing  a bright 
target.^  An  improved  predictive  track  procedure  is  described  and  illu- 
stratedvfchat  obtains  a least  squares  error  solution  to  an  atmospheric 
phase  mod^.  A visual  assessment  of  image-  degradation  with  different 
amounts  of  toase  distorti<"  has  also  been  performed  and  is  illustrated. 
While  the  im^e  qualtiy  as.  essments  are  subjective,  there  is  a correlation 
between  image\degradation  and  rms  phase  distortion. 


UNCLASSIFIED 


SECURITY  CU ASStriCATlON  OF  THIS  PAQZfWhon  D»tm  Entered) 


TABLE  OF  CONTENTS 


SECTION  PAGE 

I INTRODUCTION  TO  PREDETECTION 

COMPENSATED  IMAGING  SYSTEMS II 

II  CLASSES  OF  WAVEFRONT  ERROR  SENSORS 17 

A.  Comparison  of  Quantum  Noise  Effects  on 

the  Hartman- Test  and  Pinliole  Interfero- 
meter Systems 18 

B.  Wavefront  Correction  with  a Pinhole 

Interferometer  for  an  EIxtended  Target 33 

C.  Com^puter  Evaluation  of  Wavefront' 

Correction  with  Image  Quality  Detectors 38 

ill  DEFORMABLE  MIRROR  DESIGN 49 

A.  Surface  Figure 52 

B.  Structural  Resonances 56 

C.  Interactuator  Coupling 65 

IV  THE  MULTIPLE  APERTURE  INTERLINKED 

BLOCK  CODED  SYSTEM 77 

A.  Introduction 77 

B.  Pairwise  Element  Angle  Calculation  for 

Four  Elements  per  Block  82 

C.  Block  Rectification  Procedure  for  Four 

Elements  per  Block  . 84 

D.  Introduction  to  the  Two  Element  per 

Block  Organization 85 

E.  Element  Angle  Calculation  for  Two 

Elements  per  Block  86 


SECTION  V deleted  from  this  document 


:SECTrOfl  -VI  deleted  from  this  document 


VU  IMPROVED  PREDICTIVE  TRACK  PROCEDURES  . , . 137 

A.  One  PKase  Day ef  '" Quick  Ddok"' 

Technique  ................  137 

B.  Multiple  Phase  Layer  Predictive 

Track- Scheme  ; . . . ; . . ; . . . 140 

S.  Uhobservahi'e  Multiple  Ptese  La^^  i . . ..  . . 148 

Um  iMAGE;QUAiLTT  Y WITH  . . 15-7 

IX  REPORT  e.ON&D^^  v.  ......  . ..  171 

REFERENCES.  17 3 

APPENDIX  . , , . . . . . j . , V . .-  175 


LIST  OF  ILLUSTRATIONS 

FIGURE  ■ PAGE 

I- l  An  imaging  COAT  system  is  basically  a 

servo-mechanism 13 

II- 1 Pinhole  interferometer . 20 

II-  2 Four-quadrant  detectors  for  each  subaperture  ....  21 

II-  3 Flow  "chart  for  predetection  compensation 39 

IL-i  Phaser  diagram  for  the  electric  fields, 

using  element  by  element  method 40 

II- 5 Phaser  diagram  for  the  electric  fields, 

using  elements  by  halves  method 42 

II- 6  Walsh  functions  (1,  2,  3,  4) 44 

III-  1(a)  Sixty-four  element  I- COAT  deformable  mirror 

prototype  . . . 50 

III- 1(b)  Interferogram  of  surface  showing  one  actuator 

energized 50 

III-  2 Hughes  high  power  phase  compensating 

deformable  mirror 51 

III- 3 High  power  28  kHz  deformable  dither  mirror 

prototype 53 

TTT-4  Conceptual  phase  corrector/ dither  mirror 

configurations  . . . . ^ 54 

III-5  Faceplate -actuator  design  considerations 55 

III- 6 Isocontours  of  a portion  of  the  64- element 

I-COAT  mirror  ; . . . ' ....'.  57 

III-7  Double-pulse  holographic  interferogram  of 

the  64- element  I-COAT  mirror 58 

III-8  /.ctuator,  faceplate,  and  backpiate  resonant  modes.  . . 59 

III-9  Resonant  frequencies  of  the  faceplate/actuator 

mode,  and  actuator  column  modes  ........  6l 

lll-l  0 Interactuator  coupling  as  a function  of  ratio  of 

actuator  head  diameter  to  actuator  spacing  .....  66 


FIGURE 


PAGE 


ill- 1 1 

m-12 

III-  1 3 
III- 14 

III-  15 

IV-  1 

IV- 2 
IV-3 

IV-4 
IV- 5 


Seven  actuators  . 
P’aceplate  section 
Actuator  head 
Spring  model  . 


Linked,  block- coded  I- COAT  system . . . . . 

Typical  mirror  deformations  and  block- 

plane  images  with  piston  mirrors 


Phasor  diagram  which  defines- «}> 


67 

67 

69 

12 


Conceptual  deformable  mirror  cross-section 

using  a honeycomb  backplate  structure 

Block-coded,  multiple-apertured,  inter-linked 
imaging  COAT  system 


Linked,  block- coded  I- COAT  system 


74 

79 

79 

80 
80 
81 


FIGURES  FOR  SECTION  V AND  VI  deleted  from  this  document. 


VII- 1 Optical  geometry  for  measurement  of  phase 

layers  by  12  detectors 141 

VlI-2  Optical  setup  for  least  squares  estimation 

procedure 142 

VII-3  Phase  measurement  matrix  for  Fig,  VI.- 2 and 

3 detectors  and  9 measurement  time  intervals 145 


FIGURE 


PAGE 


VII- 4  Pseudoinverse  of  matrix  A in  Fig.  III-3 147 

VII- 5 Unmeasurable  basis  vectors  for  two  layer 

atmosphere 151 

VIII- 1 Comparison  of  computer  calculated  and 

theoretical  phase  structure  function  results 160 

VIII- 2 Photograph  of  object  on  display 162 

VIII- 3 Selective  gain  curve  applied  in  ^ABGC 

within  a window  area*  162 

VIII-4  Degraded  image,  using  phase  screen  No.  3 

for  conditions  in  Table  VIII-1 165 

VIII- 5 Degraded  image,  using  phase  screen  No.  32 

for  conditions  in  Table  VIII-1 166 

VIII- 6 Degraded  image,  smoothed  using  five 

separate  representations  of  random  phase 

screen  No.  3 for  conditions  in  Table  VIII-1  ......  167 

VIII-  7 Degraded  image,  smoothed  using  five 

separate  representations  of  random  phase 
screen  No.  32  for  conditions  in  Table  VIII-1 


168 


PREFACE 


This  final  report  was  prepared  by  Hughes  Research  Laboratories, 
Malibu,  California  under  Contract  No.  F30602-74~C-0026.  It  describes 
work  performed  from  16  September  1974  to  15  August  1975.  The  Program 
Manager  is  C.  J.  Swigert  who  assumed  this  role  in  the  last  half  of  the 
program  from  J.  A.  Jenney.  In  addition  to  the  authors  we  wish  to 
acknowledge  the  contributions  of  V.  Evtuhov,  C.  R.  Giuliano,  J.  A.  Jenney, 
L.  Miller  and  V.  Wang,  who  served  as  technical  advisors  for  the  whole 
range  of  technical  problems  that  arose  during  the  program. 


BLAMC.NOT 


I 


m- 


IN'T'RODUCTION  to  Pl^EI)F:TErTION  COMPENSATED 
IMAGING  SYSTEMS 


ii 


It  was  observed  by  the  earliest  astronomers  that  the  resolution 
achievable  with  a large  telescope  is  determined  by  atmospheric  conditions 
rather  than  by  the  size  of  the  telescope.  Generally,  it  is  observed  that  the 
resolution  limit  set  by  free-space  diffraction  theory  is  attained  only  for 
telescopes  smaller  than  a,bout  10  cm.  Although  larger  telescopes  have 
greater  light-gathering  capabilities,  they  do  not  have  Vjetter  resolution  capa- 
bilities. The  departure  from  the  predictions  of  free-space  diffraction  theory 
is  caused  by  the  fact  thet  the  atmosphere  is  not  a homogeneous  medium,  but 
rather  contains  random  density  variations  caused  by  turbulence.  The 
refractive  inde.x  variations  associated  with  these  density  variations  cause 
the  phase  and  amplitude  of  an  image  field  to  become  randomized.  Generally, 
the  perturbing  effect  of  the  phase  variations  is  dominant  and  amplitude 
scintillations  are  of  secondary  importance. 

The  objectives  of  this  predetection  compens  ited  imaging  systems 
study  are  to  determine  effective  techniques  both  for  the  measurement  of 
the  wavefront  errors  at  a telescope  aperture  and  for  the  correct. on  of 
these  wavefront  distortions.  These  techniques  are  to  be  identified  through 
analysis,  simulation,  tradeoff  studies  and  experimental  measurements  of 
the  system  elements  required  for  wavefront  error  measurement  and 
correction. 

One  approach  to  compensate  these  wavefront  distortions  is  to 
introduce  the  appropriate  (real-time)  wavefront  correction  ahead  of  the 
image-formation  process.  We  call  systems  that  achieve  such  corrections 
Imaging  COAT  (Coherent  Optical  Adaptive  Techniques)  systems.  In 
principle,  both  propagation  and  telescope  figure-induced  wavefront  errors 
can  be  substantially  reduced  by  such  systems.  Some  classes  of  COAT 
systems  may  correct  for  amplitude  variations  as  weD,  but  these  were  not 
studied  under  this  contract,  since  they  were  not  jv-lged  essential  for  long 
range  imaging. 
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Two  general  classes  of  wavefornt  error  sensors  have  been 
identified  and  explored  in  this  study.  We  have  also  examined  several  wave- 
front  perturbation  techniques  (muliiulther  techniques)  that  allow  a detector 
to  operate  in  a time  shared  mode  to  collect  the  wavefront  error  data  that 
would  normally  require  the  use  of  several  detectors  using  a non-perturbation 
mode  of  operation.  This  technique  has  been  developed  and  refined  into  the 
Multiple  Aperture  Interlinked  (MAIL)  predetection  compensation  system. 

This  technique  has  been  evaluated  in  detail  in  the  First  and  Second  Standard 
Noise  Calculations  performed  under  this  contract  and  described  below. 

Basically,  an  imaging  COAT  system,  as  illustrated  in  Fig.  I-l,  con- 
sists of  three  components:  (1)  a reflective  or  transmissive  device  which 
changes  path  length,  index,  or  phase  shift  as  a function  of  spatial  position 
over  the  beam  according  to  prescribed  error  signal  inputs,  usually  electri- 
cal signals.  Since  most  such  devices  have  a discrete  phase  control  which  is 
monotonic  over  some  subapertures  or  control  regions  called  the  elements 
(which  correspond  roughly  to  the  "elements"  of  a microwave  phase  array), 
we  call  such  a phase-correcting  system  a phaser  matrix.  (2)  An  error- 
sensing system  which  either  detects  the  departure  (in  radians)  from  uni- 
formity of  the  received  wavefront,  or  some  secondary  measure  of  its 
departure  — for  example,  local  wavefront  tilt  error.  (3)  A servo  system 
which  derives  some  measure  of  the  local  wavefront  error  and  computes  or 
establishes  the  required  error  signals  to  be  supplied  to  the  phaser  matrix. 

Early  in  the  program  we  were  motivated  to  examine  image  quality 
error  sensing  systems,  in  which  the  "wavefront"  error  sensor  operates  at  a 
system  image  plane,  preferably  the  main  image  plane  of  the  system. 

Clearly,  one  does  not  obtain  a direct  wavefront  error  measure  from  such 
information.  Alternatively,  one  must  extract  a measure  of  image  quality 
or  sharpness  based  on  on"  or  more  measures  of  image  sharpness  via  the 
image  plct-ie  detectors;  this  measure  is  then  maximized  with  respect  to  the 
available  degrees  of  freedom,  the  mirror  actuator  or  electrode  excitations. 
There  were  two  motivations  for  concentrating  on  this  approach:  (1)  the 
number  of  additional  detectors  required  for  the  approach  was  low,  perhaps 
zero;  and  (2)  confidence  in  performance  would  be  high  since  local  optical 
path  difference  errors  in  the  separate  imaging  and  wavefront  processing 
paths  following  the  beam  splitter  of  Fig.  I-l)  could  be  avoided. 
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Fig.  I-l.  An  imaging  COAT  system  is  basically  a servo-mechanism 

which  is  capable  of  sensing  phase  errors  in  the  arriving  v/ave- 
front  from  atmospheric  propagation  effect  or  telescope  figure 
distortions,  and  driving  these  to  zero  or  an  acceptably  small 
value. 


We  had  analyzed  a number  of  such  systems  on  IR&D  (prior  to  the 
present  contract)  and  for  the  sake  of  completeness  are  including  some  of 
these  results  in  the  current  report  (Section  II).  This  work  continued  only 
for  a very  brief  time  on  the  present  contract.  We  quickly  found  that  most 
such  systems  achieve  a signal-to-noise  ratio  which  typically  varies  as  the 
inverse  of  the  square  of  the  number  of  actuators  and  in  the  best  cases  as 
the  inverse  of  the  first  power.  This  remained  true  for  all  image  plane 
detection  and/or  dither  policies  of  which  we  could  conceive.  At  the  same 
time,  analysis  of  requirements  was  suggesting  that  the  number  of  actuators 
would  be  between  two  hundred  and  four  hundred  and  that  the  utmost  in 
sensitivities  would  be  required  for  applications  of  extreme  interest.  Further, 
analysis  of  the  required  speed  of  response  was  suggesting  that  it  would  be 
impossible  to  achieve  the  required  dither  excursion  speeds,  for  image 
quality  sensing,  with  conventional  PZT  driven  mirrors.  Clearly,  im.age 
plane  detection  of  the  types  considered  to  that  point  would  not  suffice.  * 

Consequently,  a new  class  of  ditn*;r  systems  was  conceived  which  we 
call  "block-encoded"  systems,  which  a oided  the  inverse  dependence  upon 
the  number  of  actuators  and  improved  detection  sensitivities  in  other  respects. 
Further,  the  required  number  of  dither  excursions  was  reduced  by  many 
orders  of  magnitude.  A particular  subclass  of  the  block-encoded  systems, 
which  we  call  Multiple- Aperture-interLinked  MAIL  appeared  to  offer  the 
best  performance  and  the  remainder  of  the  program  was  devoted  to  achieving 
the  best  performance  with  this  approach.  This  system  is  however  more 
complex  than  its  predecessors.  For  example,  a separate  detector  plane  is 
required  to  extract  the  error  information. 

Section  II  discusses  the  classes  of  wavefront  error  sensors  analyzed 
in  this  study.  Two  general  classes  of  wavefront  error  sensors,  namely 
interferometric  sensors  and  image  quality  sensors,  are  identified  and 
explored  in  this  study.  The  Searing  interferometer  and  the  pinhole  interfero- 
meter, which  is  analyzed  in  some  detail,  are  both  examples  of  interferometric 
sensors.  Image  quality  sensors  include  the  Hartman  test  system,  the  inte- 
gral of  intensity  through  a mask,  e.  g.  , a pinhole,  and  tha  integral  of  intensity 
squared  detector. 

These  systems  types  are  for  long-range  imaging  through  a turbulent  atmos- 
phere. We  feel  that  the  techniques  may  be  perfectly  suited  to  short-range 
problems  or  space  telescopes  where  long  servo  integration  times  are 
feasible. 
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Section  III  reviews  our  deformable  mirror  design  and  development 
effort  for  image  correction  and  dither.  Because  of  program  redirection  to 
emphasize  theoretical  evaluation  of  the  candidate  predetection  compensated 
imaging  concepts,  the  deformable  mirror  design  effort  was  not  carried  to 
completion. 

Section  IV  presents  the  description  of  Hughes  Research  Laboratories 
predetection  compensated  imaging  concept.  The  concept  can  be  described  as 
a block  coded,  multiple  aperture  interlinked  (MAIL),  deferred  correction, 
polystep  dith»»**  system. 

Section  V reviews  our  analysis  of  the  MAIL  compensated  imaging 
system  performance  for  the  first  standard  noise  calculation.  The  problem 
statements  for  the  first  and  second  standard  noise  calculations  were  pre- 
sented to  aid  evaluation  of  the  candidate  predetection  compensated  imaging 
system.  The  parameters  of  the  first  Standard  Noise  problem  are  sum- 
marized in  Table  V-1. 

Section  VI  reviews  our  analysis  of  the  second  standard  noise  problem 
for  the  MAIL  compensated  imaging  system.  The  phase  error  measurement 
strategy  for  the  IvLAIL  concept  is  to  use  predictive  track  of  multiple  phase 
layers  to  obtain  an  extended  integration  time,  when  it  is  advantageous; 
otherwise,  rely  on  "quick  looks"  of  the  phase  error  while  viewing  a bright 
target.  Due  to  lack  of  time,  the  predictive  track  algorithm  has  not  been 
optimally  formulated  for  the  second  Standard  Noise  calculation. 

Section  VII  presents  an  improved  predictive  track  procedure  that 
obtains  a least  squares  error  solution  to  an  atmospheric  phase  model.  The 
atmospheric  model  contains  only  two  phase  layers,  but  permits  a complete 
demonstration  of  the  track  procedure.  It  is  show  i that  there  will  always  be 
six  unobservable  phase  patterns  out  of  the  set  of  all  possible  phase  patterns 
that  can  never  be  measured  with  this  predictive  track  procedure.  Further, 
it  is  concluded  that  with  the  addition  of  more  phase  layers  to  the  atmospheric 
model,  there  will  be  more  unmeasurable  phase  patterns  added  to  the  null 
space.  However,  we  have  not  determined  the  frequency  with  which  these 
phase  patterns  will  occur. 
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degradatio,-'  with  different  amounts  of  phase  distortion.  While  the  image 
quality  assessments  are  subjective,  there  is  a correlation  between  image 
degradation  and  rms  phase  distortion. 

Section  IX  presents  the  summary  conclusions  of  the  report. 


II.  CLASSES  OF  WAVEFRONT  ERROR  SENSORS 


Two  general  classes  of  wavefront  error  sensors  have  been  identified 
and  explored  in  this  study.  The  two  generic  classes  of  wavefront  error 
sensors  are  interferometric  sensors  and  image  quality  sensors.  The 
Shearing  interferometer  and  the  pinhole  interferometer,  which  has  been 
analyzed  at  Hughes  Research  Laboratories,  are  both  examples  of  inter- 
ferometric sensors.  Image  quality  sensors  include  the  Hartman  test  system, 
the  integral  of  intensity  through  a mask,  e.  g.  , a pinhole,  and  the  integral  of 
intensity  squared  detector. 

Early  in  the  study  the  Hartman  test  system  showed  good  signal-to- 
noise  characteristics  as  a wavefront  error  sensor.  Similarly,  the  pinhole 
interferometer,  which  is  basically  M Michaelson  interferometers  in  parallel, 
with  one  for  each  of  the  M phase  shifter  elements,  has  good  signal-to-noise 
characteristics.  A comparative  signal-to-noise  ratio  analysis  of  these  two 
sensors  for  point  sources  is  presented  in  Section  II-A.  Further  analysis  of 
the  pinhole  interferometer  has  demonstrated  an  inability  to  measure  the 
correct  phase  distortion  when  viewing  multiple  point  sources  and  extended 
targets.  The  degradation,  i.  e. , increased  noise,  is  caused  by  the  mixing 
of  the  multiple  wavefronts  from  the  multiple  sources  at  each  of  the  detectors. 
An  analysis  of  this  multiple  target  process  is  presented  in  Section  II-B. 

In  parallel  with  the  analysis  of  the  Hartman  and  pinhole  interferometer, 
considerable  effort  has  been  expended  to  develop  viable  image  quality  tech- 
niques. An  analytic  computer  simulation  of  the  predetection  compensation 
process  was  developed  to  evaluate  the  performance  characteristics  of  image 
quality  improvement  measuring  (1)  integral  of  intensity  through  a pinhole 
matched  to  the  Airy  disk  size  and  (2)  integral  of  intensity  squared.  A 
description  of  the  simulation  and  the  results  are  presented  in  Section  II-C. 

The  results  indicate  the  required  target  brightness  for  good  wavefront  cor- 

2 

rection,  e.g.  , <(A$)  > less  than  0.4  or  K/10  correction,  under  typical 
seeing  conditions,  is  greater  than  the  brightness  of  most  targets  of  interest. 
An  examination  of  this  excessively  large  target  brightness  requirement 
indicates  that  the  source  of  the  difficulty  is  an  excessively  large  noise  com- 
ponent in  the  signal-to-noise  ratio  calculation.  Since  the  image  quality 
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measurement  is  made  with  only  one  detector,  the  signal  component  for  each 
element  of  the  primary  aperture  acts  like  a noise  term,  for  the  phase  mea- 
surement of  each  of  the  other  elements  comprising  the  primary  aperture. 

Our  technique  to  improve  the  signal-to-roise  ratio  is  to  use  more  detectors 
and  have  each  detector  view  only  two  or  four  actuator  elements  at  the  pri- 
mary aperture.  This  technique  has  been  developed  into  a Multiple  Aperture 
InterLinked  (MAIL)  system,  descrioed  in  Section  IV,  and  has  been  our  base- 
line concept  for  comparative  evaluation  in  the  first  and  second  Standard  Noise 
calculation  for  predetection  compensetion  imaging  systems. 


A.  Comparison  of  Quantum  Noise  Effects  on  the  Hartman-Test  and 
Pinhole  Interferometer  Systems 

This  section  mathematically  compares  the  signal-to-noise  sensi- 
tivities of  the  Hartmann  test  and  che  pinhole -interferometer  adaptive  imaging 
systems.  The  comparison  is  done  for  four  discrete  cases  and,  in  each  case, 
for  a continuum  of  targets.  The  four  cases  are 


Case  I 
Case  II 
Case  III 
Case  IV 


Filled  disk  of  n subapertures,  system  nearly 
converged 

Ring  of  n subajiertures,  system  nearly 
con\ erged 

Filled  disk  of  n subapertures,  system 
unconverged 

Ring  of  n subapertures,  system  unconverged. 


The  unconverged  state  is  that  where  the  phases  of  the  wavefrunl  in 
each  subaperture  are  nearly  stoichastically  independent,  i.  e.  , 


E 

k--i 


% 


= vn 


(2.  1) 


unconverged 


This  condition  holds,  for  example,  when  the  system  is  first  turned  on.  The 
nearly  converged  state  has  the  subaperture  phases  differing  from  '^ach  other 
by  only  fractions  of  a redian,  giving 
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(2.  2) 


k=l 

It  is  desirable  to  obtain  this  state,  most  of  the  time,  during  operation  of  the 
imaging  system. 

A continuum  of  gaussian  targets  is  considered;  the  targets  are 
assumed  to  have  a gaussian  spatial  distribution  of  reflectance  with 
standard  deviation  cr: 


= \pa 

nearly 

converged 


at  target 


-(x^+y^)/(2o-^o-^R^) 
I e \ o / 


(2.  3) 


where  R is  the  distance  to  the  target.  The  parameter  o-^  is  selected  so  that 
O'  = 1 corresponds  to  a target  just  resolvable  by  the  telescope: 


O'  = 0.  8 \/d 


(2.4) 


Table  II- 1 summarizes  the  results  described  in  this  section  for  the 
standard  deviations  o-(4)=!')  of  the  phase-error  estimators  for  the  Ha:  tmann 
test  system  and  the  pinhole  interferometer,  illustrated  in  Fig.  II- 1.  All 
results  given  in  Table  II- 1 must  be  multiplied  by  where  1.^,  is  the 

light  flux  (expressed  as  photoelectrons)  on  a subaperture. 

1.  The  Hartmann  Test  System 

The  Hartmann  test  system  has  four-quadrant  detectors  for 
each  subaperture  (Fig.  II- 2). 

The  variable  proportional  to  x-direction  tilt  is  taken  as 


I"  . I2  4I4  - I,  - I3 


(2.5) 
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Fig.  II- 1 . Pinhole  interferometer. 


20 


Table  II- 1.  Standard  Deviations  cr(4)=;'-)  of  the_J’hase-Error  Estimators 
(Multiply  all  quantities  by  l/y~Ij  where  I-j-  is  the  flux 
on  a subaperture) 


Hartmann 


Small  T 
Large  O' 
Pinhole 


Small  O' 
Large  a 


Fille 

d Disk 

Ring 

Converged 

Unconverged 

Converged 

Unconverged 

Case  I 

Case  III 

Case  II 

Case  IV 

4.  5 

4.  5 <j/v^n 

4.  5 
4.  5 O'  / 

1 . 8 \(n 
S.l  of 

1.  8 sfn 
S.l  ^ ! \/~n 

1 

0.  45  cr 

0.  45  ^/n 
0.  64  O' 

1 

0.  14  o-n 

0.  2 n 
0.6  c 

T1839 


1 2 
3 4 


Fig.  II- 2,  Four -quadrant  detectors 
for  each  subaperture. 


and  for  y-direction  tilt, 


= 'F  *2  - b - >4 


(2.6) 


(All  currents  are  in  units  of  photoelectror.s  t-^r  second.  ) The  quantum- 
noise  response  of  th"  svstem  will  be  de>'ernu.ied  by  first  determining  the 
system's  best  estirr''t  >r  for  the  phase  error  in  the  ith  subaperture,  and 
then  calculating  the  standard  deviation  of  (|>r  caused  by  the  quantum  noise. 
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Tlie  currents  h are  then  random  variables  having  roughly  gaussian 
distribution  (the  gaussian  approximation  to  the  Poisson  for  large  numbers) 
with  expectation  E(I.)  and  standard  deviation  <^(1.)  r:  ^£(1.)  » 

a . Filled- Aperture  Cases  (1  and  III) 

(1)  The  Estimator  — The  phase  error  estimates  in 
the  fllled-aperture  case  are  chosen  to  minimize  the  X*'  function 


- y^[(4>r  . -4>f  , . - .) 


(4> - kjy  ) 1 

V1,J  1,J-1  j,j/  J 


(2.7) 


where  the  subscript  i,j  refers  to  the  subaperture  in  the  i column  and  the 
jth  row  of  the  filled  aperture.  The  constant  k is  inversely  proportional  to 
the  rate  of  change  of 


witt>  change  in  <().  a small  change  A<|).  , causes  a displacem'. nt  of  the  image 
J J 

on  detector  i,  j by 


X^A4>. 


Ax  = 


(2.8) 


where  i is  the  system  focal  length  and  w is  the  spacing  between  adjacent 
subapertures.  At  the  detector,  the  image  is  appj  oximately  gaussian  with 
standard  deviation 


^ = 2cr 

image  1 


g(T  Vl/v^  + 0-^  = '/iTv^+ff 

o d 


(2.  ■)) 





Here  v is  the  ratio  of  the  full -aperture  diffraction  limit  to  that  of  a 
subaperture,  d is  the  diameter  of  the  full  aperture.  For  the  filled  disk, 

V = l//^  . {2.  10) 


The  fluxes  on  detectors  2 and  4 together  for  our  gaussian  image  are 


•CO  2 . ? / 7 7 

(u-x)^/2(f<7,)^  f (u-xr/z(u  r 
I,  + / e ^ du  / 1 e ^ du  (2.  11) 


where  x is  the  x-coordinate  of  the  center  of  the  image.  For  small  x the 
numerator  is 


I e 


“ (u-x)^/2(£o-j)^ 


du  = Ij,  I e 

/-X 


” v^/2(f  cr  j 


dv 


since 


•w  2/,,.  ,2 

V /2(fo-J 

= xl^+l^|  e dv  (2.12) 


0X 


•®  .2 
V /2(£a-^) 


'-X 


dv 


= -e 


v^/2(£o-p^ 


x=0 


= -1  . (2.13) 


So,  for  small  x, 


^2  ^4  " ~ 7^  ^ T 


ftr^  y/z^ 


(2.  14) 
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and  by  similar  reasoning  we  will  have 


Then 


" S = 


-Irj,  X 


/ /Tf 


..  2 


J 

X 


'2  "'4 


2lj  X 


(2.  15) 


The  displacement  ./2trw  of  an  initially  centered  image  then  gives  out- 

i»  J 

put  change 


AT^  - \2A4>  _ ^T 

./tt'  2itw  " T~3^r^ 

1 


^cr,  y Ztt  cr^  w(2ir)' 


(2.  16) 


The  constant  k then  is 


Al’'  ' 


3/2 


(2.  17) 


For  the  filled  aperture,  w d/yn7  so  eqs.  (2.9)  and  (2.  17)  combine  to  give 


6.  3 


2^2 
(T  + O' 

e 


cases  I and  III  ~ 


n I„ 


i).  2 Jn  + 

y«  4 


(2.  18) 


using  eq.  (2.  10)  with  cr^ 


1/v  --  yiT" 
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In  eq.  (2.  7)  X is  minimized  when  the  <|)*  ^ are  chosen  so  that 

<j>*  . = - kl^,,  . 

V i-l,J  1+1, J 1+1, J 

+ <j>.  . , + kl7  . + 4).  + kiy 

i,j  ^i,j  + l i,j+l/ 


(2.  19) 


(with  the  absence  of  appropriate  terms  at  the  boundaries). 

This  implicit  set  of  equations  can  be  solved  for  each  <)>!''  . explicitly 
• .XV  1^  1 

as  a linear  function  of  the  I 's  and  I^  's.  The  linear  operator  M giving  the 

solution  is  local  and  position- independent  (except  near  the  boundaries): 


^j-2,i-k 

k,)2 


Table  II- 2 gives  our  computer  calculation  of  M. 


(2.  20) 


Table  n-2. 

Values  of 

^ for  eq.  (2.  20) 

I ^ 

o 

' 

o 

1 

-0.  03 

-0,  02 

-0.  01 

0.  01 

0.  02 

0.  03 

0.  04 

3 

1 

o 

• 

o 

-0.  04 

-0.04 

-0.  03 

0.  03 

0.  04 

0.  04 

0.  04 

2 

-0,  05 

-0.  05 

-0.  06 

-0.  07 

0.  07 

0.  06 

0.  05 

0.  05 

1 

-0.  05 

-0.  07 

-0.  11 

-0.  25 

0.  25 

0.  11 

0.  07 

0.  05 

0 

1 

o 

o 

-0.  05 

-0.  06 

I 

o 

o 

0.  07 

0.  06 

0.  05 

0.  05 

-1 

O 

O 

1 

-0.  04 

1 

O 

o 

0.  03 

0.  03 

0.  04 

0.  04 

0.  04 

-2 

1 

o 

o 

-0.  03 

-C.  02 

-0.01 

0.  01 

0.  02 

0.  03 

0.  04 

-3 

-3 

-2 

-1 

0 

1 

2 

3 

4 

K 
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(2)  Standard  Deviation  of  — The  fluctuations  of  the 
estimators  due  to  quantum  noise  will  nov/  be  determined.  Assuming  the 
atmospheric  phase  errors  are  frozen  and  the  system  collects  photoelec- 
trons per  subaperture,  the  standard  deviations  of  the  system's  estimates 
()>!'  j of  the  atmospheric  errors  are  calculated. 

The  standard  deviations  of  each  and  are  all  \^1,^  since  each 
subape’-tare  generates  I photoelectrons.  The  standard  deviation  of  4>?'  . 
is  then  [using  Eq.  (2.20)] 


o-(4>?  .)  ^ k 7 (m^  , . . + ) 

k,£ 


k,£ 


= k [2  X 0.  26] 


since  ^ Mj^  ^ = 0.  26.  This  result  is  independent  of  the  degree  of  con- 
vergence of  the  system. 


b.  Ring  Aperture  Cases  (II  and  IV)  — The  X to  be 


minimized  is 


Here  i indexes  position  around  the  circle  and  1.  is  the  ti‘;i  in  the  tangential 

^ V 2 ^ 

direction  (compare  with  I“,  H).  This  X is  minimized  by  the  solution 


'i ' 


(2.  23) 
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where  the  kernel  Q is 


Ok" 


1 _ ii 

2 " n 


1 . ^ 1 

2 n 


l^k>0 


- 2<k^0 


(2,24) 


The  quantity  Qj^  is 


T.^‘f-2 


(2.  25) 


For  the  ring  aperture,  we  have  v = ir/n  and  w s;  ird/n.  Equation  (2.  17)  then 


gives 


, 0.  8Tr(2Tr)^/2  Vn^/tr^  + cr^  19.7  / Z,  2 , 2 

k = 5 7 = — 5 — dn  /tr  + or 

2n  n ^ 


The  standard  deviation  of  <j>‘,  using  the  same  reasoning  as  in  the 
previous  section,  is  then 


cr(cj)f)  ^ k^  /y^,Qf  " 5.  7 Vn/TT^  + g-^/n  l^fX^  , (2.26) 


again  independent  of  the  degree  of  convergence  of  the  system. 

2.  The  Pinhole  Interferometer  System 

Figure  II- 1 illustrates  the  pinhole  interferometer.  The  nar- 
row pinhole  has  diameter  d^  = 2.  44A.£/d,  the  diameter  of  the  telescope's 
Airy  disk.  Focal  length  of  the  telescope  is  f and  d is  the  diameter  of  the 
telescope  aperture.  The  wide  pinhole  has  diameter  d^  = vd^  where 
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(Z.  il) 


Jn  for  the  filled  disk 

V 

II /it  for  tiic  ring. 


We  let  6 be  the  fraction  of  light  in  each  subaperture  that  is  sent  to 
the  signal  beam,  i.  e.  , 6 is  the  reflectivity  of  Sj.  Assuming  does  not 
absorb,  the  light  in  the  reference  beam  in  the  disk  aperture  cases  is  then 
n(l  - 6)  We  will  not  treat  the  case  that  the  light  inside  the  ring  is  used 

in  the  reference  beam. 

The  amount  of  light  penetrating  the  narrow  pinhole  is  nl^.  This  is 
the  light  incident  on  the  pinhole  multiplied  by  the  fraction  of  the  gaussian 
image  falling  on  the  circular  opening  of  pinhole.  This  fraction  is  obtained 
by  a geometric-optics  calculation:  For  our  four  cases  this  is 


n I 


n 


(1  - 6)  n I.J, 

1.2  + 0.8  0-^ 


(1  - e)  nl^ 

1.  2 + 0.  8 0-^ 

(1  - 6)  n 
0.  8 + 0,  8 n 

(1  - 6)  n 

0.  8 (T^  + 0.  8 n^/iT^ 


Case  I 


Case  II 


Case  III 


Case  IV 


(2.  28) 


(The  cases  here  and  in  the  following  are  those  defined  at  the  beginning  of 
this  section.  ) 

The  amount  of  light  penetrating  the  wide  pinhole  is,  again  from  the 
geometry. 
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n I - 


w 


ne 


1 + 0.  8 0-  /n 


ne 


1 + 0.  8 cr^  iT^/n^ 


Cases  I and  III 


Cases  II  and  IV 


(2.  29) 


We  have  used  approximating  fractions  in  place  of  the  exact  functions  for  I 


n 


and  I . 
w 


a.  Depth  of  Modulations  — The  degree  of  interferometer 

modulation  I is  defined  here  as  the  peak  value  of  the  sinusoidal  variation 
m ^ 

of  interferometer  output  with  path  length  difference.  The  expected  outputs  1 
and  2 of  the  interferometer  for  channel  i are 


E 


r ’)  = I " 'w  « - V2)]  (2-  30) 


where  A is  pathlength  difference  and  v is  the  frequency  of  light.  We  will 
introduce  a fixed  it/2  path  length  difference  in  one  arm  so  that  the  zero- 
phase-error  condition  will  give  equal  outputs  in  the  two  photomultipliers. 
So  the  actual  phase  error  will  be  4>  = 2ir  v^C  Ait/2. 

The  appendix  derives  for  our  four  cases  and  for  small  tr  and 
large  cr.  Table  II-3  summarizes  these  results.  The  transition  from 
small  0”  to  large  c takes  place  between  o"  = 1 and  tr  = v. 


Table  II- 3.  Evaluation  of  I 


m 


Small  (T 


'ases  I and  II 


ll  I 
•y  w n 


Cases  III  and  IV 


i j I I 
\ w n 


Large  <r 





>'-N 


Small  (T 


Here  we  assume  cr  < 1.  Then 


<r(4'p  = 


' 0.  8 + 0.  2 e 

3.  2 e (1  - e)  I, 


0.  8 + 0.  2 e 
3.  2 e (1  - e)  Ir 


/ 1 + 0.  8 ne 

'4  e (1  - 6)  I. 


1 4 i(0.  8/it^)  - 1)  e 


46  (1  - 0)  I, 


Case  I 


Case  n 


Case  III 


Case  IV 


(2.  34) 


The  following  gives  the  0.  For  Case  III  we  alsc  give  a practical  6 assuming 
that  mirror  reflectivity  cannot  change  during  convergence.  The  resulting 
are  given. 


Small  g 

e_ 


Case  II  0.  5 

Case  III  (optimum) 

y/n 

0.  5 (practical) 
Case  IV  0.5 


o'(4>  1) 


0,  45  /iT/  J I^ 

0.  6 yfn/  J I^ 

Jo.  5(1  + 0.08  n^)/ JT^  (2.35) 


Case  I 0.  5 
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Large  (r;(r  > n 


(T  (<!>'■. ) 


1 + n e 


5.  e (1  - 6)  L 


/c/  2 2.  ^ 1 

ff  6 {n  --nr)  + iT 

2it  V 1.2  e (1  - e) 


rr^  (e  + + ne 

2.4  0(1  - 0) 


Case  I 


Case  II 


{ 


Case  III 


(T^  [g  + (1  - 0)  6/it‘ 

2.46(1  - 6)  U 


Case  IV 


As  in  the  previous  subsection,  we  give  optimum  values  of  9 for 


). 


Large  <r 


e 


o-(4>T) 


Case  I 


0.  45  cr 


Case  II 


ir/n  (no  use  of 

central  light) 


0.  14  <r  n/ 


Case  III  1 / y/~n~ 


0.64CT  /jj;^ 

0.  6 0-  n/i  + 3.  l /n/fTr  ( 


Case  IV 


■tr/n 


B.  Wavefront  Correction  with  a Pinhole  Interferometer  for  an 

Extended  Target 

Our  analysis  of  the  pinhole  interferometer  has  led  us  to  conclude  that 
there  is  insufficient  information  obtained  from  the  measurements  to  perform 
wavefront  correction,  when  viewing  an  extended  object  or  a multiple  point 
source  target.  An  assessment  of  the  analysis  given  below  is  that  too  much 
information  must  be  loiown  about  the  target  or  obtained  from  the  data  to  do 
wavefront  correction.  The  number  of  point  sources,  their  relative  bright- 
ness, and  their  positions  with  respect  to  each  other  over  the  integration 
period  must  be  known  in  order  to  estimate  the  wavefront  distortion.  These 
parameters  will  all  be  changing  as  the  satellite  rotates  in  its  orbit  and  can 
directly  introduce  significant  error  in  the  wavefront  phase  error  estimate. 

Figure  II- 1 illustrates  the  optical  setup  and  coordinates  system  for 
the  analysis  of  the  multiple  point  source  target.  Assume  the  target  is  char- 
acterized as 


N 

£(x,  y)  = 6(x  - x^,  y - y^)  (2.  38) 

n=l 

where  A^  is  the  amplitude  and  (x^,  y^)  is  the  position  of  point  source  n. 

The  amplitude  at  the  telescope  aperture  prior  to  traversing  the  phase  screen 
is  given  by 


N 


F(u,  v) 


= ^[f(x,y)]  - ^ 


A e 
n 


-if  (u- X ) 
— — n 


(2.39) 


n=l 


where  u = (u,  v),  x = (x  , y ),  f is  a scale  factor  and  Sf  is  the  Fourier 
— ' » '»  n ■'n 

transform  operator.  The  phase  screen  at  the  aperture  introduces  the  phase 
function 
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Fj(u,  V)  = F(u,  v)e 


ii{i{u,  v) 


= [f('X,  -y)  © ^(-x,  -y)] 

= [£(x,  y)  © ^(x,y)]  . (2.40) 


where  we  define  the  function 


^(x,  y) 


jr-i  [^iWu.v)] 


(2.41) 


and  note  the  Fourier  transform  of  the  product  of  two  Fourier  transform  func- 
tions is  the  convolution  of  the  original  two  functions.  The  optical  wave- 
fronts  from  the  beam  splitter  impinging  on  each  of  the  pinholes  in  planes 
and  is  given  by 


f;(x',y’)  = a.JTj-F^(u,v)] 


1 = 


1 = 


aT(x,  y)  © ^(x,  y) 


N 

a.  / A 'S'(x'  - X , y*  - y ) 
1 / j n n ^ 

n=  1 


1 for  plane  P^ 

2 for  plane  P^ 


0 < a.  < 1 
1 


(2.42) 


The  wavefront  transmitted  through  the  wide  pinhole  (assuming  no  spatial 
filtering  of  the  target  image)  to  plane  P^  is 


Fj  (u',  v')  = ^[fj  (x',  y')] 

+i2(u'  * X jj)+i‘J^(2i') 

n=l 


(2.43) 
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The  wavefroi:t  transmitted  by  the  small  pinhole  at  plane  p2»  is 


f2(x',y')  = a^  £2  (x',  y')  6{x',  y') 


^2^3 


N 


A - X , y'  - y ) 

■ n n ^ ^n 


Ln=l 


f>(x'.  y') 


The  optical  wavefront  (x')  propagated  to  plane  P,  is 

J 


F^(u',  V')  = ^[f'2{x',  y')] 


N 


= ^4  22^N*'-V  -^n' 


n=l 


N 


EB  with  B = a.  A ^(-x  , -y  ) 

n n 4 n n ^n 


n=  1 


The  optical  intensity  at  a detector  in  plane  P^  with  position  u'  = (u' 

Il(ii')  = [f;  (U-)  + F2(u-)j 


N r 

S[ 


a, A + B 
In  n 


^ +i£u'  • X +ivji 

+ a,  A B*  e " 

Inn 


-i£u 


T a, A ’"B  e 


'u'  • X -idj  1 


(2.44 


(2.45 


, v')  is 


(2.  46 


since  each  of  the  point  sources  is  incoherent  and  has  zero  mutual  intensity. 
If  we  re-express 


■ie 


'3'(-x  ) ^ e 

' n'  n 


n 


with  >1'  = 

n n 


then 


n=l 


2 ^ 2^2 
a,  + a . SE' 

1 4 n 


+ 2a, a.  cos(vb(u')  +6  + Hn' • x.  ) 

1 4 n ^ — n — 


(2.47) 


By  a similar  analysis,  the  optical  intensity  at  a detector  in  plane 
at  position  u',  due  to  a phase  shift  of  •n/2  radians  in  the  function  (u')  is 
given  by 


(u')  + F^iu')  e^^^^ 


] [f^(u')  + F'^in')  "> 


N 


= + a^4'^  + 2a^a^4'j^  cos(4j(u;)  + 6^^  - tt/2  + 2u'  • x^)j 


n=l 


N 

= R + *4  '*'R  "“1  Vn  K + 'i' 

n=  1 


(2.48) 
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If  we  design  the  optical  system  so  that  a^  = a^, 

N ^ 

I,(u')  = af  > A r [l  + + 2^  sin(4/(u')  + 0 + £u'  • x^)l 

2 — 1 / j n|L  n n ^ — n — J 

n=  1 
N 

I,(u')  = a^  A I |l  f + 2^'  cos(4'(u')+6  + fu'  * x )1 
1 — ' 1 / j n L n n 't'v_  / ^ — — n'J 

n=  1 

(2.49) 

It  would  greatly  aid  us  in  the  estimation  of  iHu')  if  we  knew  the  number  of 
point  sources,  N,  their  amplitudes,  A^,  and  their  positions  x^  = (x^,  y^), 
n = 1,  2,  . . . , N.  An  additional  difficulty  with  the  ambiguity  in  the  esti- 
mate of  4^{u')  can  be  demonctrated  easily  for  the  case  N = 1.  Assume  that 
= 1 andxj  = = 0.  Then 


l2(u')  = 

^ 2 
2a  j 

^ [1  + sin 

(2. 

50) 

« 

Ij(u')  = 

:>  2 

^ [1  + cos 

4;(u')] 

(2. 

51) 

and 


^2  , _ 1 + sin  4<(u') 

I^  — ' 1 + cos  4>(u') 


(2.  52) 


has  values  for  4^  in  the  interval  -ti/Z  < 4^  < it  and  tr  < ijj  - 3 tt/2  for  any  given 
value  of  12/1^*  The  ambiguities  become  more  complicated  as  N increases. 
Confronted  with  the  intractable  form  of  the  parameter  4“  to  be  estimated  in 
the  data  |I^,  ^2^— Hu’  have  discarded  the  pinhole  interferometer  as  a 

candidate  for  predetection  compensation. 
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C.  Computer  Evaluation  of  Wavefront  Correction  With  Image 

Quality  Detectors 

1 . Introduction 

An  analytic  computer  simulation  has  been  used  to 
evaluate  the  performance  characteristics  of  image  quality  improvement 
measuring  integral  of  intensity  through  a pinhole  matched  to  the  Airy  disk 
size  and  integral  of  intensity  squared.  The  computer  simulation  models  a 
48  in.  telescope  viewing  a point  source  (in  the  far  field)  through  a phase 
screen  directly  in  front  of  the  telescope.  Using  the  design  philosophy  of 
minimizing  system  costs  where  possible,  32  piston  corrector  mirror  ele- 
ments, arranged  in  an  annulus  on  the  primary,  relay  the  incident  illumina- 
tion to  the  single  image  quality  detector. 

Each  of  the  mirror  elements  is  10  cm  x 10  cm  in  size  and  can  be 
shifted  forward  or  backward  individually  or  in  synchrony  with  all  other  mir- 
ror elements  to  phase  up  the  distorted  wavefront  incident  on  the  primary 
aperture.  The  phase  screen  model  used  to  distort  the  wavefront  is  gen- 
erated by  a fast  Fourier  transform  technique  which  yields  a two-dimensional 
distribution  of  gaussian  random  number s having  correlation.properties  con- 
sistent with  those  associated  with  Kolmogorov  turbulence.  The  strength  of 
the  phase  screen  is  an  input  parameter  to  the  computer  simulation. 

Other  parameter  inputs  to  the  computer  simulation  include 

• The  number  of  target  generated  photoelectrons  per 
mirror  element  per  image  measurement 

• The  size  and  position  of  the  detector  in  the  focal  plane 

• Th-;  image  quality  criterion,  either  integral  of  intensity 
th-ough  the  pinhole  (if  index  SEXP  = 1)  or  integral  of 
the  squared  intensity  (if  index  SEXP  = 2). 

• An  index  parameter  specifying  that  phase  compensation 
is  to  be  performed  one  element  at  a time  in  sequence 
or  with  half  the  elements  simultaneously,  in  sequence. 

A flowchart  of  the  computer  simulation  is  shown  in  Fig.  II- 3. 
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Initialize,  read  input  parameters 


Generate  phase  screen  and  compute 
std.  deviation  of  phase  variation 


Compute  field  distribution  EL(x,y)  at 
detector  for  each  element  L=  1.....32 


Determine  E field  magnitude  of  all 
aperture  elements  over  detector  and 
plot 


Element  by  element  Choose  method  of  Half  of  elements 

] element  update,  LI=  1 ► 1 ♦ 


For  aperture  element  LI,  dither  element 
by  kir/2,  k=  -1,0,+1,  and  compute  in- 
tensity Ik(x,y)  over  detector.  Add 
poisson  noise  to  Ik(x,y) 


For  Walsh  function  W{L,LI),  L-  1,32, 
dither  element  if  W(L,LI)>0  by  kir/2, 
k=  -1,0,+1  and  compute  intensity 
Ik(>!>y)  over  detector.  Add  poisson 
noise  to  I|^(x,y) 


Compute  image  quality  criterion  for 
k=  -1,0,+1 

S(k)  = jTl(x,y)  for  Sexp  = 1 

= /l^(x,y)  for  Sexp  = 2 


Compute  image  quality  criterion  for 
k=  -1,0,+1 

S(k)  =/l(x,y)  for  Sexp  = 1 

=/l^(x,y)  for  Sexp  = 2 


Obtain  element  phase  angle  from 
S(-l),  S(0),  S(1) 

4ilt  =tan'H!sO)-s{-i))/(2s(o)-s(-u-s(i)l 


Obtain  elements  composite  phase  angle 
from 

S(-l),  S(0),  S(l) 

0|_I=  tan‘h<sn)-s'-i)/(2S(o)-s{-i)-sn)] 


Compensate  mirror  LI  by  -4>, , 
LI  = LI+1 


Have  all  elements  been  compensated^ 


. Should  we  perform  another  iteratiorC 


Compensate  each  element  L by  -<ti| , 
where  u(L,LI)>0 


Have  elements  been  fully  updated^ 


'Should  we  perform  another  iteration^ 


Are  there  more  cases  to  be  run 


Fig.  II-3.  Flow  chart  for  predetection  compensation  with  32  piston 
elements  and  image  quality  detectors. 
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2. 


Principle  of  Operation 


The  following  discussion  will  illustrate  the  principle  of  phase 
compensation  for  the  element-by-element  method  and  the  elements  by  halves 
method  for  the  integral  of  intensity  through  a pinhole  image  quality  detector. 
A discussion  of  the  simulation  results  will  conclude  the  section. 

The  element-by-element  method  of  phase  compensation  for  N ele- 
ments holds  N-1  elements  constant  and  dithers  one  element  +X-/4  and  -X./4 
wavelengths  (or  +it/2  radians)  to  estimate  the  phase  deviation  4>  between  the 
average  phase  over  the  N-1  elements  and  the  phase  of  the  remaining 
element. 

Taking  as  the  resultant  electric  field  over  the  N-1  elements  and 
as  the  field  over  the  remaining  element,  and  (f)  as  the  phase  angle 
between  the  fields  (Fig.  II-4),  the  undithered  intensity  is  given  by 

= E^  + E^  + 2EjE,j.cos4>  . (2.53) 

Shifting  the  element  forward  to  introduce  tt/2  radius  phase  shift, 

= E^  + Ej  + 2EjE.j,  cos  (ct>  + tt/2) 

= E^  + Ej  - 2EjE^  sin  4>  (2.54) 


Fig.  II-4.  Phaser  diagram  for  the  electric  fields 
using  the  element  by  element  method. 
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. 1' 
I - 


E>- 


Shifting  the  element  back  for  -ir/2  radians  phase  shift  yields 


I_j  = + Ej  + 2EjE^  sin.cj) 


(2.  55) 


Note  that  a good  estimator  for  4>  is 


$ 


tan 


-1 


'-1  - ^1 


21 


-’-1, 


(2.  56) 


The  elements  by  halves  method  of  compensation  for  N elements 
employs  as  many  measurements  as  the  element-by-element  method  but 
obtains  an  improved  signal-to-noise  ratio  by  increasing  the  modulation  of 
the  dithered  intensity.  The  method  of  implementation  chosen  here  is  to 
sequence  through  each  of  the  Walsh  functions  W^j,  LI  = 2,  . . . , 32,  for  as 
many  elements  as  are  to  be  compensated.  All  of  the  L lens  elements, 

L = 1,  . . . , 32,  will  be  shifted  +it/2  for  one  measurement  and  -tt/2  for  the 


next  measurement  if  the  value  of  the  Walsh  function  Wj^j(L)  > 0.  Other - 


A 


wise,  the  element  is  not  shifted. 

The  measured  phase  deviation  described  below,  is  used  to  com- 

pensate each  of  the  elements  L whose  Walsh  function  value  Wj^j(L)  > t. 

To  illustrate  the  by  halves  estimation  process,  assume  each  element 
has  a constant  electric  field  strength  E^  and  phasor  6^,  for  elements 
£ = 1,  2,  ...  32  (Fig.  II-5).  The  intensity  seen  at  the  detector  is 


I = 


N 


£=1 


E 


N 


N 


N + 2 


EE 


cos  (6j  - e^) 


£=1  m=£+l 


(2.57) 
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Fig.  II-5,  Phaser  diagram  for  the  electric  fields 
using  the  elements  by  halves  method. 


To  compute  the  intensity  for  step  LI  in  the  sequence,  LI  = 1,  ...»  31,  some 
elements  are  dithered  by  ir/2.  Let  be  an  index  set  with  iff  element 

is  dithered.  Let  be  an  index  set  with  iff  element  j2  i® 

dithered.  We  are  using  the  connection  "iff"  for  "if  and  only  if.  " 


K+tt/Z)  = 


N + 


- 2 


H Yj  -%) 

jj.Ji  kj.j,  V ■/ 


J €j  k €j  \ / 


'2  “2  ‘'2  2 


L E 

h'h  ^2*^2  V / 


(2.58) 


1 
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Illustration  of  the  elements  by  halves  method  is  done  most  easily 
for  the  linearized  case  with  6^  small,  j = 1,  . . . , 32.  It  can  easily  be 
shown  that  for  this  case 


I,  (+tt/2)  = N 


E 


N 

2 


I]  «jl  ■ 

Jl^Jl  j2^J2 


(2.59) 


and 


I (-it/2)  = N 


E 


(2.  60) 


Note  that  the  phase  estimator  at  small  phase  angles  yields 


N, 

4>t  = T 


jlj  (-tt/2)  - Ij  (+u/2jj  N 

jij^(-7r/2)  - Ij  (+u/2)j 

|ljVir/2)  + Ij  (+it/2)  '^^I 

X (0)  - I_  (-11/2)  - I^  (v/2) 
1 '^1  ^1 

/ N 

Jl^Jl  l2^J2 


(2. 6i; 


An  example  will  suffice  to  illustrate  how  these  measurements  are 
used  for  phase  compensations.  Consider  the  case  for  N = 4.  The  Walsh 
functions  for  N = 4 are  illustrated  in  Fig.  II-6.  For  W(m,  2),  Jj  = |3,  4} 
and  is  denoted  W(m,  3),  = \l,  3[  and  $j^  is  denoted 

For  W(m,  4),  J ^ = |l , 3{  and  $j  ^ is  denoted  We  will  also  define  a 
dummy  variable  to  ease  manipulation  of  the  variables.  For  let 
Jj  = |l,  2,  3,  4|  as  would  be  the  case  for  W{m,  1). 
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Thus,  we  obtain  ^ 


A 

4>i 

A 

4>2 

$3 

I 

-e-> 

1 

1 

1 

1 

1 

-1 

-1 

1 

1 

-1 

1 

1 

■■1 

1 

-1 

1 

-1 

CD 

1 

1 

02 

03 

1 

O 


(2.  62) 


- 1 T . T - 

By  the  orthogonality  of  Walsh  functions,  A = A , i.  e.  , 6 = A 


— 

CD 

1 

1-1-1  1 

1 

®2 

1 -1  1 -1 

h 

1111 

A 

■^3 

' 0 
1 

1 1-1-1 

A 

4>4 

and 


(2.63) 


We  can  reference  all  6 to  one  0^,  e.  g. , 6^.  By  matrix  manipulation,  we 
obtain  an  ejqiression  that  does  not  use  the  first  Walsh  function,  namely. 


■e; 

0 2 - 2 

1 

II 

2 2 0 

h 

2 0-2 

1 

(2.  64) 


Thus,  the  measured  values  $4  permit  us  to  estimate  the 

phase  angles  0.,  i = 2,  ...,  4,  for  correction.  Consequently,  given  N 
elements,  only  N-1  Walsh  functions  need  be  used  in  the  elements  by 
halves  method. 
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3.  Computer  Simulation  Results 


Tables  II-4  and  5 summarize  the  results  obtained  for  the 
cases  of  interest  that  were  examined.  Table  II-4  presents  the  results  using 
the  integral  of  the  intensity  through  a pinhole  matched  to  the  Airy  disk  as 
the  image  quality  sensor.  The  results  are  tabulated  as  a function  of  NPE, 
the  number  of  photoelectrons  generated  per  element  area  per  measurement. 
The  results  in  each  table  list  the  ratio  of  peak  intensity  on  the  detector  to 
the  maximum  possible  intensity,  prior  to  correction  and  after  one  complete 
iteration  to  update  all  32  elements. 

The  trends  of  the  results  are  evident  for  both  the  element-by-element 
method  and  the  elements  by  halves  method.  The  element-by-element  method 
does  not  appear  to  converge  (with  Poisson  noise  at  the  detector)  for 
NPE  < 500  photoelectrons /element  area/measurement  with  either  image 
quality  sensor.  The  elements  by  halves  method  appears  to  work  signifi- 
cantly better  for  both  image  quality  sensors,  being  able  to  converge  until 
NPE  is  reduced  to  100  or  150  photoelectrons /element  area/measurement. 
These  results  are  not  cause  for  optimism  for  the  elements  by  halves  method 
since  it  does  require  a large  number  of  photoelectrons /iteration  to  update 
the  elements. 


Table  II-4,  Ratio  of  Peak/Maximum  Possible  Intensity  With  No 
Correction  and  After  One  Iteration  Using  the 
Integral  of  Intensity  Squared  Detector.  Standard 
deviation  of  phase  screen  = 1.  17  radians.  NPE  is 
number  of  photoelectrons /1 00  cm^  element/ 
measurement. 


NPE 

Element  by  Element  Method 

N Elements  by  Halves  Method 

1000 

0.  34 0.  50 

0.  34 0.  62 

500 

0.  34 -0.  33 

0.  34 0.  52 

250 

— 

0.  34 0.50 

125 

— 

0.  34  — ► 0.44 

100 

— 

0,  34 '0.  29 

75 

— 

0.  34 0.  22 

1 


Table  II-5.  Ratio  of  Peak /Maxi mum  Possible  Intensity  Prior  to  Correction  and  After  On< 
Iteration  Using  the  Integral  of  Intensity  Through  a Pinhole  Detector.  NPE  is 
number  of  photoelectron/ 1 00  cm^  element/measur ement. 


As  a brief  example,  the  minimum  number  of  photoelectrons  for  one 

iteration  of  32  elements  with  the  elements  by  halves  method  can  be  taken  as 

(31  Walsh  functions)  (3  measurements /Walsh  function)  (100  PE/element  area/ 

measurement)  - 9300  photoelectroncs/iteration  with  the  elements  having  an 

2 

element  area  of  100  cm  . This  required  number  of  photoelectrons  for  this 
method  is  much  higher  than  the  number  available  for  many  of  the  objects  of 
interest  and  is  cause  for  us  to  attempt  to  reduce  this  required  photoelectron 
count.  Our  improved  method  of  aperture  correction,  multiple  aperture 
interlinked  (MAIL),  is  described  in  Section  IV. 


III. 


DEFORMABLE  MIRROR  DESIGN 


In  i.ae  initial  period  of  this  compensated  imaging  program,  a portion 
our  effort  was  placed  on  the  design  of  deformable  mirrors  for  image 
correction  and  dither.  This  effort  was  not  carried  to  completion  because 
of  program  redirection  to  emphasize  theoretical  evaluation  of  the  candidate 
predetection  compensated  imaging  concept.  Deformable  mirror  design  is 
emphasized  because  it  offers  batter  phase  error  correction  for  a given 
number  of  actuators  than  the  segmented  mirror  and  it  also  allows  the 
incorporation  of  a number  of  performance  features  not  compatible  with 
monolithic  designs.  The  degree  of  development  of  deformable  mirror  tech- 
nology makes  it  the  first  choice  for  phase  correction  in  our  system  concepts. 
Because  of  its  superior  voltage  sensitivity,  accurate  approximation  to  the 
desired  surface  figure,  compatibility  with  simple  cooling  arrangements  (if 
required  by  actuator  dissipation),  and  state  of  development,  the  deformable 
mirror  driven  by  a set  of  discrete  actuators  is  our  prime  candidate  for  a 
combined  dither /phase  corrector  at  present. 

The  following  topics  will  be  discussed  in  this  section: 

• Surface  figure 

• Resonant  modes 

• Interactuator  coupling 

• Assembly  techniques 

Figure  III- 1 is  a S4-element  ejqjerimental  Imaging-COAT  prototype 
mirror  assembled  prior  to  this  contract  to  test  several  aspects  of  a simple, 
high  actuator  density  design.  A thin,  relatively  flexible  glass  faceplate  is 
supported  by  an  array  of  actuator  heads,  which  are  in  turn  supported  by 
hollow  cylindrical  PZT  actuators.  Figure  III- 1(b)  shows  the  mirror  figure 
as  obtained  after  initial  polishing.  Several  polishing  techniques  were 
tested,  but  no  attempt  was  made  to  make  this  a thermally  stable  structure. 

Figure  III- 2 shows  a deformable  mirror  delivered  to  the  Navy  by 
Hughes  for  phase  correction  at  10.6  pm  in  a high  power  laser  system.  For 
high  power,  cooled,  deformable  mirrors  for  10.6  pm,  the  continuous 
surface  mirror  actuated  by  discrete  piezoelectric  or  hydraulic  actuators 
is  the  favored  approach. 
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Figure  III-3  is  a high  power  dither  mirror  developed  on  a parallel 
IR&D  program.  In  this,  nine  actuator  stations  are  included,  with  the  center 
driven.  The  faceplate  is  grooved  to  relax  faceplate  stiffness.  An  interesting 
feature  of  this  design  is  its  extremely  high  operating  frequency  (over  40  kHz), 
made  possible  by  a cylinder  spring  surrounding  the  PZT  cylinder  to  raise 
its  resonant  frequency,  and  the  symmetrical  model  attachment,  which 
decouples  the  dither  drive  oscillations  from  the  backplate,  giving  better 
isolation  between  actuators. 

The  three  basic  types  of  phase  correctors  considered  for  compen- 
sated imaging  applications  are  schematically  illustrated  in  Fig.  III-4.  The 
segmented  array  of  Fig.  III-4(b)  features  certain  simplicities  of  analysis, 
and  its  performance  readily  lends  itself  to  an  optimized  dither  mirror,  as 
will  be  discussed.  The  monolithic  array  (Fig.  III-4(c))  has  strong  advan- 
tages in  economy  of  manufacture  and  maximum  array  density.  Balanced 
against  this  is  a relatively  undeveloped  technology  and  distinct  disadvan- 
tages in  high  voltage  drive  requirements  and  the  lack  of  cooling  provisions. 
The  third  type  of  phase  corrector  (Fig.  III-4(a))  is  the  deformable  mirror, 
discrete  actuator  corrector.  We  discuss  the  design  considerations  for 
the  deformable  mirror  below. 

A.  Surface  Figure 

The  manner  in  which  a surface  approximates  the  desired  figure  by  a 
finite  number  of  actuators  significantly  influences  the  system  performance. 

It  appears  that  some  considerable  advantage  can  be  gained  by  using  a surface 
which  can  optimally  tilt  as  well  as  displace. 

The  deformable  faceplate  configuration  might  be  expected  to  provide 
an  inherent  tilt  correction.  However,  to  achieve  the  optimum  correction, 
the  actuator  faceplate  interface  must  be  carefully  configured  to  optimize  the 
mirror's  ability  to  fit  a typical  atmospheric  dislocation  function,  while 
maintaining  faceplate  stress  within  acceptable  bounds.  A few  sample 
configurations  are  shown  in  Fig.  III-5.  Figure  IIl-5(a)  shows  the  nearly 
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REFLECTING  SURFACE 


RIBBON  CABLE  TO  AMPLIF'ERS 


c Monolithic  mirror 

Fig.  III-4.  Conceptual  phase  cor  rector/dither  mirror  configurations 


segmented  behavior  of  a deformable  mirror  when  actuator  head  diameters 
(large  compared  to  the  spacing)  are  used.  The  opposite  extreme  (Fig,  III-5(b)) 
can  result  in  high  localized  stress  and  distortion  over  the  head.  Figure 
III-5(c)  proposes  a conceptual  actuator  head  which  provides  flexure  to  accom- 
modate tilt,  with  minimum  compromise  in  inter -actuator  coupli;i7.  Fig- 
ure IIl-5(d)  illustrates  a simpler  approach  where  a relatively  thick  faceplate 
is  first  fabricated,  and  then  portions  are  etched  away  to  increase  faceplate 
compliance.  This  approach  is  promising  because  it  allows  careful  shaping 
of  the  faceplate  cross  section  to  achieve  the  desired  optimum  mirror  figure. 

In  our  work  to  date  we  have  investigated  deformable  mirrors  both 
analytically  and  experimentally.  The  effect  of  deformable  mirror  figure 
upon  COAT  performance  has  been  modeled  in  computer  simulations  where 
the  faceplate  deflection  has  been  approximated  as  the  deflection  of  a circular 
disk,  clamped  at  the  edges,  and  deflected  at  a central  point.  An  analytic 
solution  to  this  model  is  expressed  as  a table  of  deflections  which  is  then 
incorporated  into  a COAT  simulation  algorithm. 

The  central  analytical  approximations  to  the  surface  figure  assumed 
a circularly  symmetric  smooth  displacement  from  an  actuator.  Experi- 
mental observation,  however,  revealed  a more  complex  pattern.  Figure  III-6 
shows  a static  interferometric  measurement  of  faceplate  displacement  for  a 
hexagonal  array  of  actuators  on  a deformable  mirror,  A double-pulse 
holographic  exposure  of  the  same  mirror  being  driven  at  1 3 kHz  is  shown  in 
Fig.  Ill- 7,  giving  the  identical  pattern  as  for  the  static  excitation. 

We  have  shown  that  the  actual  surface,  with  finite  head  diameter  and 
fixity  at  six  points  around  the  actuators,  as  described  in  Fig.  Ill- 6,  can  also 
be  modeled  with  our  computer  techniques.  This  requires  a simple  change 
of  the  table  used  in  our  computer  algorithm.  As  a result,  a more  accurate 
measure  of  deformable  mirror  COAT  performance  is  obtained. 

B.  Structural  Resonances 

The  key  structural  resonances  in  our  deformable  mirror  system  are 
the  actuator /faceplate  piston  mode  and  the  backplate  mode  (Figs.  III-S(al 
and  III-8(b)).  The  frequency  of  the  actuator /faceplate  mode  is  determined  by 
the  mass  of  the  faceplate,  the  mass  of  the  actuator,  and  the  spring  constant 
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Fig.  III-7.  Double-pulse  holographic  inter- 
ferogram  of  the  64-elemem  I- 
COAT  mirror  with  one  actuator 
driven  at  13  kHz,  ±250  V. 
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RESONANT  MODES 
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FACEPLATE  DRUMHEAD  MODES 


Fig.  III-8.  Actuator,  faceilate,  and  backplate  resonant  modes. 


of  the  actuator.  In  this  system,  the  coupling  to  adjacent  actuators,  via 
the  faceplate,  contributes  a small  amount  to  the  overall  stiffness.  The 
effect  of  the  coupling  is  to  simply  broaden  the  actuator  resonance  into  a 
band  of  resonances. 

Additional  resonant  modes  that  must  be  considered  are  lateral 
column  resonances  (Fig.  ni-8(c))  and  faceplate  interactuator  drumhead 
modes. 

1 . Actuator /Faceplate  Resonance 

The  lowest  resonant  frequency  occurs  when  all  the  actuators 
vibrate  in  unison,  so  that  the  faceplate  moves  as  a large  piston.  This 
resonant  frequency  can  be  obtained  from  the  following  expression: 


f 


o 


J_ 

2it 


M + m/3 


(3.  1) 


where 

k = spring  constant  of  actuator 
m = mass  of  actuator 

M = effective  mass  of  faceplate  defined  as  the  mass  of  a disc 
whose  radius  is  two  thirds  of  the  actuator  spacing 

For  an  actuator  1 cm  long,  0.  317  cm  diameter,  and  0.  76  mm  wall  thick- 
ness, we  obtain 


f^  = 40  kHz  . (3.2) 

Figure  III-9  shows  the  effect  of  the  actuator  length  on  the  resonant 
frequency.  The  resonant  frequency  given  by  this  curve  will  not  be  affected 
appreciably  by  changes  in  actuator  wall  thickness  or  diameter  since  changes 
in  these  quantities  leave  the  ratio  of  k/m  constant.  In  addition,  changes  in 
the  faceplate  mass  will  have  a nearly  negligible  effect  on  the  resonant  fre- 
quency since  the  faceplate  mass  is  small  compared  to  the  actuator  mass. 
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2. 


Inter  actuator  Faceplate  Resonant  Modes 


A stiff  faceplate  will  raise  the  resonant  frequency  of  the 
unsupported  faceplate  region  between  actuators  but  a flexible  faceplate  will 
nainimize  the  necessary  actuator  forces  and  interactuator  coupling.  The 
resonant  frequency  of  this  region  is  given  by  solutions  of  the  form 


f 

n 


(3.  3) 


where 

E = Young's  modulus  of  elasticity 
p = density 


V = Poisson's  ratio 


r = radius  of  disk 


t = thickness 


For  "I  rectangular  array  on  0.4  cm  centers,  the  lowest  resonance 
can  be  approximated  by 


° K(0.‘4|2) 


(3.4) 


For  a 0.  4 mm  thick  glass  faceplate  we  get  a resonant  frequency  of 
about  125  kHz.  In  other  words,  it  appears  the  faceplate  can  be  made  as  thin 
as  optical  techniques  permit. 

3.  Lateral  Column  Resonances 

These  modes  are  orthogonal  to  the  driving  forces  as  well  as 
well  as  to  the  direction  of  phase  shift.  Thus  their  effect  upon  corrector  mir- 
ror performance  is  certainly  smaller  than  the  previous  cases,  and  perhaps 
neglectable.  Nevertheless,  an  estimate  of  these  frequencies  can  be  easily 
made,  where 
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for  the  lowest  mode,  I,  shown  in  Fig.  III-8(c). 
where 

I = moment  of  inertia  of  cylinder  cross  section 
)JL  = W/gi,  mass/unit  length 
^ = length 

f =12  kHz  for  a 1 cm  actuator,  0.  317  cm  diameter. 

° \ 

Figure  III-9  illustrates  the  ^rong  inverse  square  length  dependence. 

\ 

However,  since  the  faceplate  can  be  *iBtrongly  constrained  at  the  edges  against 
lateral  movement  (in  a deformable  bu^not  a segmented  mirror),  the  mode 
pictured  in  Fig,  III-8(c),  II,  may  be  mWe  appropriate.  For  this  mode,  the 
resonant  frequency  for  our  typical  actuator  is  then 


f s 50  kHz  . (3. 6) 

n 

This  type  II  mode  is  above  the  passband  for  actuators  up  to  2 cm. 

In  addition,  use  of  a 0.  5 -cm  diameter  actuator  will  raise  the  frequency  by  a 
factor  of  2.  5,  placing  the  lowest  mode  above  the  passband  for  2 cm  actuator 
lengths. 

4.  Backplate  Resonances 

When  the  actuator  array  size  becomes  large,  structural 
resonant  modes  of  the  backplate  can  occur  at  low  frequencies  and  can  domi- 
nate faceplate  and  actuator  modes.  The  wave  equation  for  simple  thin  plates 

v4  , ^ . 0 (3.7) 

Et  dt 
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where 

r)  = displacement 
p = density 
V = Poisson's  ratio 
t = thickness  of  plate 
E = Young's  Modulus 

can  be  solved  for  circular  plates  clamped  at  the  edge. 

for  the  lowest  mode  and  for  the  edges  free,  clamped  in  the  center 
(umbrella  mode^ 


f 


n 


0.  17 


V ) 


(3.9) 


(The  overtones  in  this  system  are  not  harmonic.  ) 

For  a 15  cm  diameter  glass  substrate,  4 cm  thick,  the  lowest  fre- 
quency obtainable  (in  the  umbrella  mode)  is 

f^  = 6700  Hz  , (3.  10) 

Thus,  one  should  employ  a mirror -mounting  system  which  will  con- 
strain the  mirror  to  a higher  normal  mode,  or  use  a structure  of  higher 
specific  stiffness.  For  instance,  by  appropriate  choice  of  mounting  points, 
it  should  be  fairly  simple  to  constrain  the  mirror  motion  to  a mode  which  is 
at  least  as  high  as  the  clamped  mode  (f^  = 18  kHz). 

Alternately,  since  f^  improves  only  by  ^E/p,  materials  choice 
offers  limited  advantage;  experience  in  optical  table  design  has  shown  that 
far  greater  results  can  be  achieved  without  compromise  of  other  properties 
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(such  as  coefficient  of  thermal  expansion)  by  the  use  of  high  stiffness 
honeycomb  construction.  For  glass,  this  has  been  sometimes  implemented 
by  fusing  an  array  of  tubes  to  two  faceplates.  Corresponding  resonant 
frequencies  can  easily  be  tripled  by  this  technique. 

C.  Interactuator  Coupling 

A property  of  deformable  mirrors  that  has  been  a cause  for  concern 
is  the  extent  to  which  the  motion  of  one  actuator  affects  the  mirror  figure  at 
neighboring  stations.  It  is  shown  below  that  even  for  an  actuator  of  4 cm 
in  length  and  0.  317  cm  in  diameter  the  coupling  coefficient  is  <1%  for  a 
0.  4 mm  faceplate.  The  effect  of  faceplate  thickness  is  summarized  in 
Fig.  Ill- 10.  This  chart  plots  interactuator  coupling  and  stress  as  a function 
of  actuator  head  diameter  and  faceplate  thickness.  The  conclusions  to  be 
reached  here  are  that  faceplate  thickness  in  the  thinned  etched  regions 
(Fig.  III-5(d))  are  to  be  limited  to  about  0.  5 mm  by  considerations  of  stress. 
If  a material  other  than  glass  (or  quartz  or  Cer-Vit)  is  to  be  used,  one  with 
a lower  modulus  of  elasticity  and/or  higher  stress  limit,  \:hicker  faceplates 
can  be  used.  For  actuators  which  are  shorter  than  the  baseline  1 cm  length, 
and  therefore  stiffer,  the  interactuator  coupling  is  proportionately  reduced. 

The  interactuator  coupling  is  assumed  to  occur  only  through 
the  faceplate.  The  coupling  is  expressed  as  a percentage  of  the  initial 
motion,  and  by  linear  superposition;  the  general  case  of  arbitrary  dis- 
placement of  all  actuators  can  be  modeled.  Frequencies  are  assun^ed 
to  be  sufficiently  below  the  lowest  faceplate  resonance,  that  static 
analysis  is  valid. 

Figure  III- 1 1 shows  seven  actuators  of  a hexagonal  array.  The 
spring  forces  on  the  central  actuator  arise  from  the  faceplate  and  can 
be  approximated  by  considering  the  faceplate  to  be  terminated  at  the 
first  ring  of  surrounding  actuators.  This  section  of  faceplate  can  be 
further  approximated  as  a circular  disc,  with  a hard  peripheral  sup- 
port and  a center  load  via  the  driving  actuator  (in  Fig.  Ill- 12).  The  dis- 
placement at  the  center  of  this  model  is 
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Fig.  111-11. 
Seven  actuators. 


W = TOTAL  APPLIED  LOAD 
E = YOUNG'S  MODULUS  = 10^  psi  FOF  GLASS 
t = thickness  OF  PLATE 

m = (ur’=  (POISSON'S  RATIOr^  =4  FOR  GLASS 
a = RADIUS  OF  DISC  = INTERACTUATOR  SPACING  = 0.4  cm 
b = RADIUS  OF  ACTUATOR  HEAD  DISC 

REF  ROARK,  RELN  NO.  22  P223 


Fig.  ni-12.  Faceplate  section. 
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g I deflections  are  much  smaller  than  the  thickness.  Thus  the  load 

I deflection  and  load-stress  relations  are  linear. 


As  a sample  case,  let  us  evaluate 
t =0.4  mm  (faceplate  thickness) 
a =0.4  cm  (actuator  spacing) 
b = 0.  1 cm  (head  diameter  = Z mm) 

then 


y = 2.^^xlO”^W  in  English  units 


or,  for  y = IISO  nm  = 4.  53  x 10~^  in. 


m 


w =0.21  lb 


The  maximum  faceplate  stress  is  given  approximately  by 


where 
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2a^(m  + 1)  log  ^ + a^(m  - 1)  = b^(m  - 1) 
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whichever  is  larger. 
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Since  the  second  expression  is  larger,  the  maximum  stress 
occurs  over  the  actuator  head  (Fig.  Ill- 13)  and  is  equal  to 

<r^  = 3, 400  psi 
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Fig.  Hi- 13.  Actuator  head. 


This  is  a value  which  can  be  accommodated  by  conventional  glasses. 

For  a faceplate  which  is  properly  etched,  the  stresses  can  be  expected 
to  be  lower  due  to  the  smooth  fillets. 

The  faceplate  spring  constant  is  simply 

= 4.  64  X 10^  Ib/in.  (3.  14) 

By  comparison,  the  actuator  spring  corstan'"  can  be  calculated 
using  the  expression 


y 


[AiEJ  A^E, 


(3.  15) 
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where 


^ = actuator  length 

Aj  = cross-sectional  area  of  actuator  tube 

Ej  ~ moduli' H of  elasticity  of  PZT  in  Y^j  direction 
when  s -t  circuited 

h = height  of  actuator  head 

= average  cross  sectional  area  of  actuator  nead 

= Young's  modulus  of  elasticity  of  actuator  head 

For  our  sample  case,  let  us  take  the  more  challenging  case  of  a long 
actuator  configured  for  low  voltage  (since  this  aids  in  driver  design), 
whe  re 

actuator  length  = 4 cm  = 1.57  in. 
actuator  tube  wall  thickness  = 0.020  in. 
actuator  diameter  = 0.125  in. 

Ej  = (12.3  X 10”^^  m^/nt)"^  for  PZT-4,  short  circuited 
- 8.  1 X 10^^  nt/m^  = 11.8  x 10^  psi 
h = 0.4  cm  = 0.  1 57  in 
Aj  = 1.  5 X 10'^  in.^ 

-2  2 

A^  = 1.068  X 10  in.  , assuming  hollow  heads  of 
0.040  in.  wall  thickness 

E^  = 10  X 10^  psi  for  glass 


y ^ 1.07  X 10”  W 


*K  - 9.  34  X 10^  Ib/in. 


This  should  represent  very  nearly  a worst  (lowest  stiffness)  case 
analysis  for  the  actuator. 


70 


The  stress  in  the  actuator  due  to  the  faceplate  constraint  is 


<T  = 


W 

A 


14.  5 psi 


which  is  well  below  the  8000  psi  limit  set  by  the  manufacturer  to  hold 
stress  induced  depoling  to  10%.  Dynamic  stresses  may  however  be 
higher. 

The  buckling  limit  for  this  actuator,  incidentally,  is 


W 


tt^EI 


where 

1/4  4\ 

I = -^  IT  I R - 1 the  moment  of  inertia 

R = O.D.  4 2 

R = I.  D.  4 2 

0 

-6  4 

1 = 9, 42  X 10  in. 

Whien 

W = 442  lb,  or  at  stress  of  30  ksi,  in  the  actuator 
(failure  is  at  75  ksi) 

The  interactuator  coupling  can  be  deduced  from  the  following 
considerations. 

1.  The  force  on  the  actuator  head  is  equal  and  oppo- 
site to  the  force  of  constraint  of  the  faceplate. 

2.  The  force  of  a central  actuator  is  exactly  encoun- 
tered by  an  opposite  force  on  the  six  surrounding 
actuators.  Thus,  the  force  on  one  adjacent 
actuator  is  1/6  the  force  on  the  deployed  actuator 
(see  Fig.  Ill- 11). 

3.  In  Fig.  Ill- 14,  a spring  model  for  Fig.  Ill- 11  is  con- 
structed, k^(Ax.)  = 6ka(Ax2),  by  statement  (1) 

and  (2),  where  Axj  is  the  displacement  of  the 
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(a) 


(b) 

Fig,  III- 14.  Spring  model. 


faceplate  and  activated  actuator,  and  Ax2  is  the 
displacement  of  the  adjacent  actuator.  It  follows 
that  the  coupling  between  two  actuators,  expressed 
as  a fraction  is  then 


A^2 

A^l 


6k 


8.4x10 


-3 


10. 


Since  computer  simulations  have  shown  that  considerable  interaction 
coupling  can  be  tolerated,  the  faceplate  can  be  made  up  to  twice  the 
thickness  used  here.  Even  with  this  increase  in  thickness  the  inter- 
actuato!'  coupling  is  less  than  7%.  However,  faceplate  stresses  will 
have  to  be  carefully  calculated  and  controlled  by  actuator  head  diameter 
choice  and  the  geometrical  shaping  of  the  juncture  Whereas,  in  this 
example,  there,  is  only  a 5%  shortening  of  the  actuator  by  the  faceplate 
constraint,  a stiffei  faceplate  will  require  an  actuator  with  an  appre- 
ciably greater  free  extension. 
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Computer  simulations  for  multidither  COAT  systems  have  shown 
that  interactuator  coupling  of  up  to  30%  results  only  in  a small  loss  of 
convergence  time.  Thus,  with  coupling  typically  under  1%,  interactuator 
coupling  alone  does  not  significantly  impact  our  baseline  design. 

1 . Assembly  Techniques  for  Deformable  Reflectors 

During  the  past  seven  years  Hughes  Aircraft  Company  has 
gained  considerable  experience  in  joining  complex  mirror  structures.  Our 
molybdenum  cooled  mirrors  for  high  power  applications  are  now  joined  by  a 
high  temperature  brazing  process,  follov/ed  by  cold  cycling  down  to  80°K, 
to  relieve  braze  stresses  and  stabilize  the  mirror  figure.  The  deformable 
reflector  sho'vn  in  Fig.  Ill- 2 employs  threaded  preloading  rods  which  act  to 
join  faceplates  to  PZT  actuator,  preload  the  PZT-5H  stacks,  and  also  to 
supply  a portion  of  the  actuator  spring  construction.  This  design  ha.s  also 
proved  to  be  stable  with  time. 

One  prototype  for  imaging  applications  (Fig,  III-l)  uses  epoxy 
(Armstrong  A- 12)  bonding  between  the  faceplate  and  actuate.*  heads,  as  well 
as  for  the  rest  of  the  structure.  The  stability  of  such  organic  bonds  is 
presently  under  investigation.  Other  adhesives  presently  under  investigation 
are  Eastman  type  910  adhesive,  a UV  curing  PMM  Loctite  glass  adhesive, 
and  various  Owens-Illinois  silicon  glass  adhesives.  Other  techniques  that 
might  merit  consideration  include  optical  contacting,  metal  diffu.  ion  bond- 
ing, and  soldering. 

One  Hughes  preliminary  design  is  illustrated  in  Fig.  Ill- 15.  A 
15  cm  backplate  around  4 cm  thick  of  Pyrex  or  Cer-Vit  honeycomb  construc- 
tion is  mounted  to  a blocklike  mirror  mount  in  a manner  which  minimizes 
resonances.  A discrete  actuator  array  of  500  PZT-4  or  PZT-8  cylinders 
between  0,  5 cm  and  2 cm  in  length  is  ground  and  lapped  to  equal  length  in 
one  batch.  These  are  bonded  to  the  optically  flat  backplate,  with  the  aid  of 
a locating  fixture  plate.  Actuator  heads  of  a material  with  identical  expan- 
sion coefficient  to  the  faceplate,  with  a head  diameter  of  about  half  the 
actuator  spacing  are  fabricated  as  shown,  bonded  to  the  PZT  tubes,  and 
ground  and  lapped  to  create  a flat  surface.  A faceplate  of  roughly  1 mm 
thickness  is  made  optically  parallel  using  a technique  which  has  been 
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Fig.  III-15.  Conceptual  deformable  mirror  cross- section  using 
a honeycomb  backplate  structure. 
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employed  by  our  optical  fabrication  shop  in  which  the  faceplate  is  bonded  to 
a flat  during  polishing.  The  bark  surface  is  then  photomaskcd  and  etched  to 
thin  the  interactuator  regions  to  0.4  to  0,5  mm.  This  plate  is  then  bonded 
to  the  actuator  heads  and  the  final  structure  is  thermally  cycled,  stabilized, 
"aged"  by  driving  the  PZT  actuators,  and  given  a final  polish  to  the  best 
possible  flatness  if  necessary.  Possible  added  features  are  cooling  passages, 
and  an  RTV  low  modulus  damping  material.  One  thousand  addressable  leads 
are  required  <n  the  baseline  design,  and  2500  in  the  digital  design.  These 
are  appropriately  routed  through  connectors  on  the  back  of  the  mirror  mount, 
and  through  ribbon  cables  to  the  electronic  driver  rack. 
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THE  MULTIPLE  APERTURE  INTERLINKED  BLOCK 
CODED  SYSTEM 


A.  Introduction 

This  section  presents  the  description  of  the  block  coded,  multiple 
aperture,  interlinked  (MAIL),  deferred-correction,  polystep-dither 
imaging  system  proposed  for  predetection  compensated  imaging  of  objects 
viewed  through  turbulent  atmosphere.  The  linked,  block-coded  system 
dithers  a mirror  surface  sequentially  in  a code  organized  by  blocks  of  the 
mirror  surface.  The  linked,  hlock-coded  system  uses  the  polystep  dither 
sequence  to  measure  the  atmospheric  phase  distortions  introduced  in  the 
incident  wavefront  from  the  target. 

The  procedure  to  estimate  these  phase  distortions  is  as  follows. 
The  light  that  is  reflected  from  each  of  the  "elements"  (mirror  regions) 
within  a block  is  focused  to  an  image  plane.  Each  block  has  an  associated 
pinhole  in  the  image  plane  that  transmits  a portion  of  the  incident  light  to 
an  optical  detector  behind  the  pinhole  (see  Fig.  IV-1).  The  elements  of 
the  mirror  surface  may  be  either  piston  mirrors  or  portions  of  a deform- 
able mirror.  The  blocks,  lenses,  pinholes,  and  detectors  are  replicated 
many  times,  such  that  the  entire  aperture  (or  its  image)  is  completely 
covered.  Each  actuator  basically  perturbs  the  illumination  fields  falling 
on  two  adjacent  lenses  so  that  displacements  in  any  one  elemental  area 
of  the  mirror  displace  (or  shape)  the  images  in  two  adjacent  block  image 
planes.  (With  deformable  mirrors  the  perturbation  effects  may  extend 
over  more  than  two  lenses,  but  these  are  basically  second  order  effects. 
As  the  dither  process  is  sequenced  through  a series  of  discrete  steps 
(typically,  four  or  six),  the  associated  intensities,  as  transmitted  through 
the  pinhole,  are  detected,  integrated  and  stored  in  a computer  memory. 
Conceptually,  at  least,  the  associated  integrated  intensities  or  photoelec- 
tron counts  are  first  employed  to  correct  the  intrablock  errors  for  every 
block.  Since  dual,  generally  conflicting  estimates  of  the  corrections  are 
obtained  for  the  full  array  of  actuators,  a rectification  routine  is  required 
to  remove  or  at  least  minimize  the  discrepancies. 


^CSDINS  PASS  BLAmt-NOT^PIIMEB 


Many  different  dither  sequences,  block  organizations,  and 
rectification  routines  are  possible.  One  of  the  most  promising  block 
organizations  is  with  four-element  blocks  employing  a six-step,  double- 
element dither  sequence,  as  illustrated  in  Fig.  IV-2  and  described  in 
Table  IV -1.  Each  of  these  dither  sets  perturbs  the  "block"  image  which 
is  formed  by  a particular  block  lens  upon  a pinhole  located  in  its  focal 
plane,  as  illustrated  (by  typical  perturbation  sets)  in  Fig.  IV -3.  In  general, 
the  perturbations  either  displace  the  image  or  distort  it,  such  that  an 
increased  or  reduced  number  of  photons  impinges  on  the  underlying 
detector.  Fig.  IV-4.  This  increased  or  decreased  photon  covint  is  caused 
by  the  mirror  displaceme  nr  iDcin  g an  improvement  or  degradation  of  the 
image.  The  atmospheric  phase  distortions  measured  by  the  variation  in 
dithered  image  intensity  is  then  compensated  by  mirror  correction. 

A second,  promising  block  organization  employs  two  element  blocks 
using  a single  two  step  dither  sequence  as  illustrated  in  Fig.  IV -5  and 
described  in  Table  IV-2.  This  block  organization  offers  an  improved 
signal-to-noise  ratio  for  viewing  dim  targets  but  requires  many  more 
detectors,  assuming  the  spacing  between  actuators  is  constant. 


Table  IV-1.  Six-Step  (Arc-Sine)  Block  Dither  Sequence  for  a 
Four-Element  Block-Coded  Linked  I-COAT  System 
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^0 

4-1 

4^2 

4^3 

1 

0 

0 

tt/2 

.it/2 

2 

0 

0 

-tt/2 

-V2 

3 

0 

tr/2 

tt/2 

0 

4 

0 

-ir/2 

-tt/2 

0 

5 

0 

t:/2 

0 

TT  /2 

6 

0 

-tt/2 

0 

- tt/2 
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Fig.  IV- 1.  Block-coded,  multiple-apertured, 

inter-linked,  imaging  COATsystetn, 


ELEMENT 


Fig.  IV-2.  Linked,  block-coded  I-COAT 
system.  Two  elements  are 
simultaneously  dithered  in  every 
block. 
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MIRROR  DEFORMATIONS  AS  SEEN  BY  A BLOCK  LENS 


MAGNIFIED  BLOCK  PLANE  IMAGES 

Fig.  IV-3.  Typical  mirror  deformations  and  block-plane 
images  with  piston  mirrors. 
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Fig.  IV-4.  Phasor  diagram  which 
defines  . 
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Fig.  IV-5,  Linked,  block  coded,  I-COAT  syste’^n.  Adjacent 
actuators  dither  in  alternate  sequence. 
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Table  IV -2.  Two-Step  (Arc -Sine)  Block  Dither  Sequence  for  a 
Two-Element  Block-Coded  Linked  I-COAT  System 


n 


0 


4> 


1 


-tt/4 

+tt/4 


+tt/4 
- Tr/4 


B.  Pairwise  Element  Angle  Calculation  for  Four  Elements 

per  Block 

Consider  the  piston-mirror  dither  policy  of  Table  IV -1  and  suppose 
that  the  transmitted  opticc*!  intensities  through  the  pinholes  are  proportional 
to  the  on-axis  ii->tensities  (this  is  certainly  true  for  small  pinholes  but  can 
also  be  shown  to  hold  for  larger  pinholes)  within  a scale  factor.  Let  the 
fields  of  elements  0,  1,  2,  and  3 be  designated  by  the  phasors  Uq,  Uj, 
and  Uj.  Let  be  the  phasor  sum  of  and  Next,  consider  the 

eleme  it  pairings  leading  to  Uq^  and  U23  for  example.  In  steps  N = 1 and 
2,  the  pair  of  elements  "2"  and  "3"  are  perturbed  as  an  entity  in  phase 
steps  +it/2  and  -ir/2,  giving  on-axi.3  intensities  at  the  pinhole. 


4 = k[| 


u 


oil 


+ 2 


^01  ^23 


cos(  + tt/2  j + 
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(4.  1.  a) 


I2  = K 
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Oil 


+ 2 


^01  ^23 


cos(d)^  - tt/2)  + 


IJ 


23l 


(4.  i.b) 


where  K is  an  optical  constant,  relatable  to  the  incident  flux,  and  is  the 
phasing  error  of  the  phasor  pair  relative  to  (see  Fig.  IV-4). 

These  field  intensities  are  transmitted  through  a pinhole  with  trans- 
mission efficiency  Hp  and  integrated  on  a photosurface  generating  a photo- 
current with  quantum  efficiency  q . These  currents  are  then  integrated  for 
the  time  period  T^,  producing  a mean  total  number  of  photoelectrons  NPE, 
where  for  the  near  convergence  condition,  each  of  the  optical  intensities  of 
(4.  1)  gives 


4 


f 
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and  (NPE)^  is  the  mean  efi»  ctive  number  of  photoelectrons  induced  by  the 
photons  falling  on  one  (square)  element  in  one  measurement  integration 

tim  e,  T ; or 

' m 


(NPE)  = 2KU^,  T = n,  n P 
e 01  m p q t 


(4.  2.  c) 


whe^e  Is’q  is  the  number  of  photoeleclrons  per  actuator  clement  in  one 
measurement  integration  time  and  P^  is  the  total  number  of  photons  falling 
on  an  element  during  an  integration  time.  Equation  (2.  2c)  is  essentially  a 
statement  of  the  conservation  of  photons  (or  energy)  in  the  imaging  process. 
Based  on  these  photoelectron  counts,  one  forms  an  estimator  on  called 

A A 


_ J (NPE)^  - (NPE)j 


. - i 
i jn 


(NPE)^  + (NPE)j 


, 2(NPE)  sin  d'  + N - N, 
- 1 e A 7 1 


2(NPE)  + N,  + N, 
e 2 1 


where  and  N,  are  zero-mean  noise  disb’.rbances  with  variances  asso- 

C J 

dated  with  the  integrated  photoelectron  measurements  of  (4,2),  For  the 
high  S/N  case  where  and  are  very  small,  (4.  3)  becomes,  to  a 
sonable  approximation. 
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Similarly,  we  estimate  two  other  "pairwise"  phase  errors  <j)g  and 
based  on  photoelectron  counts  from  dither  sequences  N = 3,  4,  5,  and  6 


‘-k. 
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•'1 
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4>B  = s-.n 


(NPE)^  - (NPE)3 


4>^  = sin 


(NPE)^  + (NPE)3 


^ (NPE)^  - (NPE)^ 
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(4.  5) 


(4.6) 


The  true  pairwise  phasing  errors  are 


= ^^23  - *01  " 2 ^*^3  + ‘♦’2  ■ *1  ■•‘^0^ 


- *02  = 2 (*3  +'^1  -‘*’2  -*0> 


= '^12  ■ *03  = i^*2  '^'^’1  ■ *3  ■ *^0^ 


(4.  7.  a) 


(4.  7.b) 


(4.  7.  c) 


v/here  d>  is  the  phase  associated  with  the  field  U . The  right-hand 
expansion  of  (4.  7)  employs  the  property  (unique  to  two  elements)  that  the 
composite  phase  (error)  of  a phasor  pair  is  the  average  of  the  individual 
phases  (errors).  Equation  (4,  7)  yields  a set  of  estimation  correction 
settings  (within  each  block). 

C.  Block  Bectification  Procedure  for  ' our  Elements  per  Block 

Using  the  results  from  the  previous  section,  the  phase  variation 
between  elements  in  a single  block  is  obtained  from  (4.  7)  to  be 


<^>1  -*0  = *B'^'^C 


‘^'2  - ‘^’O  " *A 


(4.  8.  a) 
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(4.8.b) 


‘•’'-t-'i'i  I'll' 
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(4.8.  c) 


If  the  telescope  required  only  one  block  for  wavefront  correction,  actuators 
1,  2 and  3 could  be  displaced  to  drive  4>  ^ q»  " ‘I’ q zero.  However, 
the  telescope  requires  many  contiguous  blocks  with  each  actuator  (mirror 
element)  shared  between  two  blocks.  Wavefront  correction  for  one  block 
may  conflict  with  the  mirror  element  displacement  requirement  for  the 
element's  adjacent  block.  The  method  chosen  to  resolve  these  potential 
conflicts  is  to  select  a predetection  compensation  performance  criterion 
to  be  optimized  and  then  determine  the  actuator  control  law  that  will  seek 
to  optimize  this  performance  criterion.  We  have  chosen  to  minimize  the 
sum  over  all  the  actuator  locations  of  the  squares  of  the  optical  path  length 
errors.  The  specific  description  of  this  criterion  and  the  optimum  actuator 
control  law  is  discussed  in  Section  V-B,  It  is  shown  that  the  optimum  block 
rectification  procedure  is  of  the  form 


X = B ^ 


(4.9) 


where  ^ are  the  measured  phase  variations  B is  a matrix  obtained 


from  the  noise  covariance  matrix  of  the  measurements  and  X are  the  best 


estimate  of  the  actuator  positions  {Xj^}.  The  expression  for  B is  described 
in  Section  V-B  to  obtain  a minimum  sum  of  the  squares  of  the  optical  path 
length  errors. 


Introduction  to  the  Two  Element  per  Block  Organization 


Another  promising  block  organization  uses  two-element  blocks  in 
a two-step,  two-element  dither  sequence,  as  illustrated  in  Fig.  IV-2  and 
described  in  Table  IV-2.  Each  of  the  dither  steps  perturbs  the  "block" 
image  which  is  formed  by  a particular  block  lens  upon  a pinhole  located  in 
its  focal  plane.  The  perturbations  displace  the  image  such  that  an 
increased  or  reduced  number  of  photons  impinges  on  the  underlying 
detector.  This  increased  or  decreased  photon  count  is  caused  by  the 


mirror  displacement  being  an  improvement  or  degradation  of  the  image. 
The  atmospheric  phase  distortions  measured  by  the  variation  in  dithered 
image  intensity  is  then  compensated  by  mirror  correction. 


E.  Element  Angle  Calculation  for  Two  Elements  per  Block 


Consider  the  piston-mirror  dither  policy  of  Table  IV -1  and  suppose 
that  the  transmitted  optical  intensities  through  the  pinholes  are  proportional 
to  the  on-axis  intensities  (tnis  is  certainly  true  for  small  pinholes  but  can 
also  be  shown  to  hold  for  larger  pinholes)  within  a scale  factor.  Let  the 
fields  of  elements  0 and  1 be  designated  by  the  phasers  Uq  and  U^.  In 
steps  N = 1 and  2,  the  elements  "0"  and  "1"  are  perturbed  in  phase  steps 
+it/4  and  -it/4  alternately,  for  a net  difference  of  ir/Z,  giving  on-axis 
intensities  at  the  pinhole 
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cos((|)  + tt/2)  + 


(4.  10a) 
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(4.  10b) 


where  K is  an  optical  constant,  relatable  to  the  incident  flux,  and  (fc  is  the 
phase  error  of  the  phaser  relative  to  U^. 

These  field  intensities  are  transmitted  through  a pinhole  with  trans- 


mission efficiency  1]^  and  integrated  on  a photosu-^face  generating  a photo- 


current with  quantum  efficieiicy 


These  currents  are  then  integrated  for 


the  time  period  T^,  producing  a mean  total  number  of  photoelectrons  NPE, 
where  for  the  near  convergence  condition,  each  of  the  optical  intensities  of 
(4.  10)  gives 


(NPE)j 
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(NPE)  2 

= (NPE)^ 
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(4.  11b) 
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and  (NPE)  is  the  mean  effective  number  of  photoelectrons  induced  by  the 
photons  falling  on  one  quarter  of  an  actuator  (square)  element  in  one  mea- 


surement integration  time,  T^;  or 


N, 


(NPE)  = ZKU;:  T P 

' e 0 m 'p  'q  e 


(4.  lie) 


where  Nq  is  the  number  of  photoelectrons  per  actuator  element  area  per 
measurement  integration  time  and  P^  is  the  total  number  of  photons  falling 
on  a quarter  of  an  actuator  element  during  an  integration  time.  Equation 
(4.  11c)  is  essentially  a statement  of  the  conservation  of  photons  (or  energy) 

in  the  imaging  process.  Based  on  these  photo-electron  counts,  one  can 
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form  an  estimator  on  called  4> 
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(4.  12) 


where  and  are  /ero-mear  noise  disturbances  with  variances  asso- 
ciated with  the  integrated  photoelectron  measurements  of  (4,  11).  For  the 
high  S/N  case  where  N2  and  are  very  small,  (2.  3)  becomes,  to  a 
reasonable  approximation. 


^ = 4> 


(4.  13) 


1-^ 


For  the  special  case  when  the  target  is  of  uniform  brightness  and 
target  scintillation  is  not  present,  then  only  one  dither  step  is  required  to 
measure  6.  For  the  targets  being  considered  in  this  analysis,  the  mean 


effective  number  of  photoelectrons  (NPE)^  is  constant  and  can  be  estimated 
from  prior  data.  Hence,  the  phase  error  4»  can  be  estimated  by 
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VII. 


IMPROVED  PREDICTIVE  TRACK  PROCEDURES 


A critical  technology  issue  identified  for  the  conditions  specified  in 
the  Second  Standard  Noise  Calculation  for  any  predetection  compensated 
imaging  system  is  to  compensate  a fast-moving  image  viewed  through  strong 
plus  distributed  turbulence.  The  phase  error  measurement  strategy  for 
compensation  is  to  use  predictive  track  of  multiple  phase  layers  to  obtain 
an  extended  integration  time  when  it  is  advantageous;  otherwise,  rely  on  a 
one  phase  layer  "quick  look"  of  the  phase  error  while  viewing  a bright 
target. 

A.  One  Phase  Layer  "Quick  Look"  Technique 

The  residual  phase  error  for  our  compensated  imaging  system 
comes  from  three  sources;  theerror  introduced  by  detector  noise;  the 
error  introduced  by  an  imperfect  fit  of  the  deformable  mirror  surface  to 
the  desired  correction;  and  the  error  introduced  by  the  fact  that  the  required 
phase  correction  changes  as  a function  of  time,  and  hence  a detector  having 
a finite  response  time,  tends  to  measure  a time-averaged  value  of  the 
required  correction  rather  than  the  instantaneous  value.  Since  these 
errors  are  statistically  independent,  the  net  mean  square  residual  phase 
error  is  equal  to  the  sum  of  the  individual  errors 

<(^cb)^>  = <(A<f>)^>j^  + < + <(Acj))2>^^  , (7.1) 

where  the  subscripts  D,  M,  and  TA  denote  detector,  mirror,  and  time 
average,  respectively.  Utilizing  results  described  in  reference  A,  we  can 
express  the  terms  .n  eq.  (7.  1)  as  follows 
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E 3 A,  , 5/3 
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where  A is  the  mirror  area  per  actuator,  t is  the  measurement  time,  P 
is  the  number  of  photoelectrons  detected  per  unit  area  of  the  primary  mirror 
per  millisecond,  r^  is  the  phase  coherence  length  defined  by  F ried,  and 
c(cd)  is  given  by  the  expression 


= 2. 91  k‘ 


As  system  designers,  we  wish  to  select  the  element  area  A^  and  the 
system  response  time  t so  as  '.^o  minimize  the  net  residual  phase  error.  To 
determine  the  optimum  values  of  these  parameters,  let  us  express  eq.  (7.  1) 
in  the  form 
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The  minimum  of  <(A<{))^  with  respect  to  t and  A^  occurs  when  9<(A4>)  >/9t 
and  a<(Ad>)^>/9A  are  zero.  This  yields  the  conditions 
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which  yields  the  relations 
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with 
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(7.  14) 


To  solve  these  equations,  we  first  solve  eq.  (7.  14)  numericall-y  and 
then  utilize  the  resulting  value  of  p in  eq.  (7.  12)  to  obtain 

The  following  analysis  of  the  rms  phase  error  is  presented  for  the 
conditions  of  a satellite  at  300  km  altitude,  with  an  atmospheric  trans- 
mission cf  257o,  a quantum  efficiency  of  10%  and  an  optical  transmission 
efficiency  of  33%.  The  value  of  c(®)  and  a are  obtained  using  the  phase 
structure  function  specified  in  the  Second  Standard  Noise  Calculation.  The 
parametfrrs  for  the  rms  phase  error  obtained  by  optimally  solving  for 
and  the  measurement  time  t are  presented  in  Table  VII- 1. 

However,  it  should  be  noted  that  the  rms  phase  error  can  be  improved 
by  treating  the  atmospheric  turbulence  as  a single  phase  layer  at  altitude 
h^,  translating  at  velocity 

The  new  value  of  c(°°)  becomes 


A 

Ch  (-) 

o 


2.  91  k 


dh  C^(h) 


w + 6 (h  - h ) 
o 


3/3 


(7.  15) 
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Table  VII-1,  One-Look  RMS  Phase  Error  Measurement 

Parameters,  b =0,  = 7.  06 

o 


Target 

Intensity, 

w/sr 

P, 

photo  electrons 
per  cm^ 
per  ms 

(A  ) 
e opt 

T .,  ms 
opr 

N,  II  of 
actuators  for 
D = 1 50  cm 

2000 

50.  52 

19.97 

0. 0496 

0.  482 

885 

1000 

25.  26 

26.  55 

0.  0571 

0.  543 

666 

500 

12.  63 

35.  29 

0.  0658 

0.  611 

500 

250 

6.  315 

46.  88 

0.  0755 

0.686 

377 

125 

3.  158 

62.  24 

0.  0868 

0.  771 

284 

62.5 

1.  579 

82.  62 

0. 0998 

0.  867 

1 

214 

For  the  value  of  h that  minimizes  ^ (oj),  » = 2.  66,  reduced  from  a = 7 . 06. 
o o 

The  optimized  parameters  for  the  rms  phase  error  for  these  new  conditions 
are  presented  in  Table  VII-2. 


ir 

'ig 


B.  Multiple  Phase  Layer  Predictive  Track  Scheme 


The  phase  layer  estimation  procedure  for  the  second  Standard  Noise 
Calculation  assumed  a linear  array  of  12  detectors  viewing  the  satellite 
through  12  contiguous  columns  of  atmosphere  (Fig.  VII-1).  In  the  analysis 
below,  a least  squares  multiple  phase  layer  estimation  procedure  is  pre- 
sented for  an  array  of  M detectors  viewing  the  satellite  for  N unit  time 
intervals  (Fig,  VIl-2).  A unit  time  interval  T is  the  time  required  for  a 
phase  patch  at  the  maximum  altitude  of  interest  to  move  one-half  of  a lens 
v/idth.  This  interval  T is  the  same  unit  integration  time  used  in  the  second 
Standard  Noise  Calculation.  It  is  also  assumed  that  each  of  the  phase 
measurements  at  a detector  is  obtained  in  the  unit  time  interval  T. 

The  multiple  phase  layer  procedure  described  below  can  be  used 
for  an  arbitrary  number  of  pnase  layers.  An  error  analysis  is  also 
presented  for  this  estimation  procedure.  To  simplify  the  analysis,  a 
complete  solution  is  presented  for  a two  phase  layer  atmosphere.  The 
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Table  VII-2,  One-Look  KMS  Phase  Error  Measurement 
Parametets,  a = 2.  66 


Target 

Intensity, 

Lj/sr 

P. 

photoe! actrons 
per  cm^ 
per  ms 

(Ae)  X 
e^opt, 

cm^ 

T , ms 
opt' 

N,  II  of 
actuators  for 
D = 150  cm 

2000 

50.  52 

15.  70 

0.436 

1126 

1000 

25.  26 

20.  88 

hh 

0.492 

846 

500 

12.  63 

27.  75 

0.  105 

0.  55tt 

636 

250 

6.  315 

36.  88 

0.  121 

0.  621 

479 

125 

3.  158 

43.  98 

0.  139 

0.  698 

361 

62.5 

1 

1.  579 

65.  04 

0.  159 

0.  785 

272 

error  analysis  below  identifies  six  different  types  of  phase  patterns  that 
cannot  be  measured  by  the  detector  array,  Three  of  these  phase  patterns 
will  not  degrade  the  image  and  the  other  three  patterns  probably  do  not 
occur  often. 

The  phase  model  used  for  the  non-overlap  geometry  of  Fig.  VII-2 
assumes 

A4.(i,  t)  = +|,  l)  - 4,Ji  - i,  l)  J 

[*t  ■^1'  i) 't’t(^  i) 

where  Acj)(i,  t)  is  the  phase  measurement  A<|>  of  thei  detector  ( -4, 

-2,  0,  2,  4,  . . , ) at  time  t and  4^(i  a,  p)  is  the  lumped  parametex  alue 
of  the  phase  patch  at  time  t and  position  (i  + a,  Ph^i:  h^  the  altitude  of  the 
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top  phase  layer.  It  is  also  useful  to  note,  for  formulation  of  the  measure- 
ment process  that 


4>  = 'i>j.(i+<3'  + r.p,  p ) . 


(7.  17) 


In  the  formulation  of  the  least  squares  problem,  we  are  given  a 
real  m x n matrix  A of  rank  k < min  (m,  n) 


A£  = 


(7.  18) 


obtained  from  (7.  16)  and  (7.  17)  and  the  set  of  measurements  we  chose  to 
make.  Figure  VII-3  illustrates  matrix  A for  a choice  of  19  measurements 
made  over  nine  time  intervals  at  three  detectors.  Given  any  real  m-vector 
Ac^,  we  wish  to  find  a real  n-vector  ^ that  minimizes  the  euclidean  length 


||Ai-  i||.  The  solution  to  the  least  squares  problem  obtains  the 

orthogonal  decomposition 


A = H R K 


(7.19) 


where  H is  an  m x m orthogonal  matrix,  R is  an  m x n matrix  of  the  form 


Rll  0 


0 0 


(7.  20) 


R, , is  a k X k matrix  of  rank  k,  and  K is  an  n x n orthogonal  matrix.  It  is 
a theorem  that  the  unique  minimum  length  least  squares  solution  is  given 


A 

^ = K 


Rll  0 


0 0 


^ Aj  Ai 


A^ 


(7.21] 


144 


Fig.  VII-3.  Phase  measurement  matrix  for  Fig.  VI-2  and  3 detectors  and 
9 measurement  time  intervals. 


ft 


The  pseudoinverse  Aj  for  A in  Fig.  VII- 3 is  shown  in  Fig.  VII-4. 
The  rank  of  A is  16;  A^is  a 19-vector  and  ^ is  a 22-vector,  Since  the 
rank  of  A is  six  less  than  22-vector  .{>,  there  are  six  basis  vectors,  i.  e. , 
independent  phase  patterns,  that  cainot  be  obtained  from  the  measurements 

• 

As  a check  that  matrix  Aj,  the  pseudoinverse  of  A is  correctly  com- 
puted, we  can  proiect  forward  the  values  of  (|>  to  obtain  the  estimated  values 
of  A$(i  + 2,  t + 5T),  where 

A4.(i+2,  t+5T)  . - 2J-,  l)  - - 3 IM 


(7.22) 


from  (7.  16)  and  (7.  17).  Forming  the  propagator  matrix 

= (0,  0,  0,  0,  0,  0,  0,  1,-1,  0,  0,  0,  0,  0,  0, . 5, . 5,  -.  5,  -.  5,  0,  0,  0) 

(7.23) 

A 'T 

where  A<(>(i  + 2,  t + 5T)  = P ^ and  ^ is  defined  in  Fig.  VII-3.  It  Cein  be 
shown  that 


Mi 


A|(i+2,t+5T)  = P^J  = P^Aj  A^ 


. *(i-2i,  a)  - + (i-3i,  i) 

-T'  l) 


which  is  mathematically  identical  with  definition  (?.  ';2b 
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(7.24) 
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Similar  techniques  can  be  used  to  compute  the  rms  phase  error  of 
the  estimate  for  A4>(i  + 2,  t + 5T),  the  phase  correction  at  time  t + 5T  after 
all  measurements  are  concluded.  By  computing  the  phase  correction 
estimate 


A$(i  + 2,  t + 5T)^  = Aj  A^ 


(7.2t>) 


the  variance  of  A$(i  + 2,  t + 5T)  is  given  by 

< aJ,  a4 (i  + 2,  t + 5T)>  = <{P’^  Aj  A^)'^,  p"^  Aj  A^> 

= P^  Aj(P^  Aj)^  <A$^,A$  > 

= 1.75  <(A<j>)^>  (7.26) 

for  the  three  detectors  and  19  measurements  of  Fig.  VIII-3,  The  use  of 

more  detectors  over  more  time  intervals  will  considerably  reduce  the 

A 2 

predicted  rms  phase  error  <(A<}>)  >. 

C.  Unobservable  Multiple  Phase  Layers 

All  matrices  for  the  two-layer  phase  model  that  have  been  evaluated 
in  this  analysis  have  always  been  reduced  to  a six  dimensional  subspace  of 
xinobservable  phase  patterns.  In  order  to  determine  whether  these  unobserv- 
able phase  patterns  C2in  badly  degrade  the  image  and  to  assess  whether 
these  unobservable  phase  patterns  are  fundamental  to  the  two  layer 
measurement  process  or  are  just  artifacts  of  the  phase  estimator  design 
procedure,  it  was  felt  necessary  to  identify  a set  of  six  basis  vectors  for 
this  null  space.  The  basis  vectors  are 

1.  4>  ^i ij  = a^  i = 0,  ±1,  ±2,  ...  (7.27) 

j~0,  "fc.  5,±1,±1.5,... 
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<j> 


1 if  j+.  5 integer 
0 if  j integer 


(7.34) 


Examples  of  each  of  the  six  basis  phase  patterns  are  given  below. 

The  patterns  are  given  for  a single  point  in  time  with  the  top  layer  phase 
pattern  appropriately  positioned  above  the  lower  layer  phase  pattern  appro- 
priately positioned  above  the  lower  layer  phase  pattern  (Fig.  VlI-5). 

Having  identified  the  six  basis  vectors  of  the  null  space,  it  becomes 
important  to  determine  whether  or  not  these  unobservable  phase  patterns 
are  a fundamental  element  of  the  two  layer  measurement  process.  It  is 
claimed,  on  the  basis  of  the  following  demonstrative  argument,  that  these 
six  phase  patterns  are  fundamentally  unmeasurable  using  the  M contiguous 
phase  detectors  operated  along  parallel  lines  of  sight  for  N sets  of  time 
measurements. 

The  first  pattern  |<^(i  + 1/2,  1)  = and  <}>(j  + 1/4,  1/2)  = 0 | is 
not  measurable  because  for  all  i 

<!>  (i  +^,  ij  - t (i  -I-,  l)  = 0 (7.35) 

The  second  phase  pattern  |<}>(i -t-  1/2,  1)  = 0 and  <|>(j  + 1/4,  1/2)  = b^  | 
i^  not  measurable  because  for  all  j 


The  third  phase  pattern  |<|>(i  + 1/2,  1)  = 0 and  <}>(j  + 1/4,  1/2) 
" (-1)^^  b^  J is , not  measurable  because  for  all  j 


(7.  36) 


.s  |<j»  (j  • ~)  ' '»•  (i  • 7)  - '!>(.)  -7. 1)  - (j  - |»r)  I 


= .5(-l)' 


■'  ^o{<- 


1)2  + (-1)^  - (-1)°  - (-1)"^}  = 


n.  37) 


The  fourth  phase  pattern  ](b(i  + 1/2,  1)  = 0 and  (|»(j  + 1/4,  1/2) 
■1)^^  [j]  is  not  measurable  because  for  all  j 

= . 5(-l)^^  b^  I j - j - (j  - 1)  + (j  - 1)|  = 0 for  j integer 


. 5(-l)^j  b 


for  j + 2'  integer. 


(7.  38) 


The  fifth  phase  screen  | <f>{i  + 1/2,  1)  - i • a^  and  4>(j  + 1/4,  1/2) 
[ j ] [ is  not  measurable  because  for  all  i 


4,  (i  +|,  1)  - <{»  (i  - J,  1)  = 


i*a  -{i-l)a  = a 
o 00 


(7.  39) 


and  for  all  j 


^ ^ ^o  ( [j  r]  1 ■ [j  " 2]  ■ N - 1 ] ^ = - % 


(7.40) 
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The  contributions  and  -a^  from  each  phase  layer  combine  at  the  detector 
to  provide  zero  measured  phase  difference. 

The  sixth  phase  pattern  { 4>  (i  + 1/2,  1)  = [ i/2]  a^  and 
*(jo  + j + 1/4,  1/2)  = (-1)^^‘’‘^|  [ j]  ([  j]]  - 1)  - <j>  I 1®  measurable 
because  for  all  i 

/ a for  i even 
= 1 ° 

^ 0 for  i odd  (7. 41) 

When  i is  odd,  is  an  integer  and  when  i is  even,  + . 5 is  ein  integer. 

The  lower  layer  phase  contribution 


- + j -?■  r)-  + j t)  I 

I 

I 0 for  j integer 

I (-l)^io  a^  for  j + . 5 integer 


Setting  + . 5 an  integer  when  i is  even  synchronizes  the  two  phase  layers 

so  the  contributions  a and  -a  combine  at  the  detector  on  alternate  time 

o o 

steps  to  provide  zero  measured  phase  difference.  On  adjacent  time  steps, 
each  phase  layer  contributes  zero  phase. 
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It  has  been  shown  in  this  two-phase  layer  error  analysis  that  there 
will  always  be  six  unobservable  phase  patterns  and  they  are  identified. 

With  the  addition  of  more  phase  layers  to  the  atmospheric  model,  there 
will  be  more  unmeasurable  phase  patterns  added  to  the  null  space.  If  this 
estimation  procedure  is  pursued,  those  phase  patterns  should  be  identified. 
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VIII.  IMAGE  QUALITY  WITH  RMS  PHASE  ERROR 


We  have  carried  out  computer  simulation  studies  of  the  imaging 
COAT  system  to  determine  if  any  correlation  exists  between  image  quality 
and  residual  rms  phase  error  of  the  segmented  mirror.  For  the  simula- 
tions, a 1.  5 m diameter  segmented  mirror  was  considered  which  was 
capable  of  correcting  for  both  phase  and  tilt  distortions  in  an  incident 
wave.  Two  different  segmented  mirrors  were  considered  in  the  simulation 
studies,  differing  in  the  size  of  the  mirror  segments.  That  is,  the  nximber 
of  mesh  points  encompassed  by  each  mirror  segment  was  different.  For 
the  coarsely  segmented  mirror,  each  segment  encompasses  31  mesh  points, 
while  for  the  finer  segmented  mirror  14  mesh  points  were  included  in  each 
segment.  A brief  description  of  the  phase  and  tilt  correction  code  of  the 
segmented  mirror  is  as  follows.  The  phase  and  tilt  of  each  segment  was 
adjusted  to  give  a least  squares  fit  to  the  phase  of  the  distorted  field  over 
that  segment.  The  difference  between  the  least  squares  fit  and  actual 
phase  distribution  was  introduced  as  a residual  phase  error  in  the  system, 
i.  e.  , the  aperture  phase  distribution  in  the  PSF  determination  was  equated 
to  this  residual  phase.  Comparison  with  results  of  more  complicated 
simulations  in  which  deformable  mirrors  are  driven  by  a COAT  servo  syetem 
indicates  that  the  least  squares  technique  utilized  in  these  studies  yields 
comparable  results.  ^ 

The  phase  screen  used  in  the  simulation  studies  is  generated  by  a 
fast  Fourier  transform  technique  which  yields  a two-dimensional  distribution 
of  gaussian  random  numbers  having  correlation  properties  consistent  with 
those  associated  with  Kolmogorov  turbulence.  In  particular,  we  represent 
the  phase  in  the  form 

^(x)  = fdK  a(K)  exp  (i  K • x)  (8.1) 


where  4>,  is  the  spatial  spectrum  of  the  phase  variation,  K is  a spatial  wave- 
<P  _ “ . 

number,  and  a(K)  is  a completely  random  function  with  zero  mean  and  unity 

variance,  i.  e. , 
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<a>  = 0 


<a(Kj)  a'‘'(K^)>  = OKj  ^ K^  , 


= = ^^2 


Note  that  with  this  representation,  the  random  function  4 (^c)  satisfies  the 
WeinerKhinchin  theorem 

<4»(xj)  4.*(x2)>  = fdK  $^(K)  exp[i  K ♦ (Xj  - x^)]  , (8.^ 


as  is  required  for  any  valid  representation  of  a homogeneous  random 
1 

process. 


for  $, 


in  these  studies  we  have  utilized  a modified  von  Karman  spectrum 


0.490 

$ (K)  = 2 9 , , 

~ r 2 


where  K is  2it  divided  by  the  outer  scale  of  turbulence  L , and  r is  the 
o o'  o 

phase  coherence  length  defined  by  Fried.  The  phase  structure  function 
associated  with  this  spectrum  is 


D^(r)  = 


,-5/3  [3  <«o  r)  ] 

^o)  . 5 f(TT75) J 


where  is  a modified  Bessel  function  of  the  5/6th  order.  In  the  limit 

of  an  infinite  outer  scale,  the  phase  structure  function  given  in  eq,  (8.6) 


reduces  to 


D^(r)  = 6. 88  (r/r^)' 
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which  is  the  phase  structure  function  specified  in  the  standard  noise 
calculation  problem  statement. 

To  verify  that  our  numerically  generated  phase  screens  do  indeed 
satisfy  the  theoretical  relation  given  in  eq.  (8.  6),  we  have  calculated  the 
phase  structure  function  associated  with  an  ensemble  of  100  realizations  of 
phase  screen  sdamples  generated  by  a fast  Fourier  transform  evaluation 
of  eq.  (8.  1),  A 64  X 64  mesh  was  used  in  these  calculations,  and  it  was 
assumed  that  themesh  points  are  1.  106  cm  apart  and  that  the  outer  scale  of 
turbulence  is  2 m.  The  results  of  the  calculation  are  shown  in  Fig.  VIII- 1 
along  with  a plot  of  the  theoretical  expression  given  in  eq.  (8.  6)  and  the 
corresponding  infinite  outer  scale  expression  given  in  eq.  (8.  7). 

The  results  given  in  Fig.  VIII- 1 yield  the  following  observations: 

(11  Our  numerical  technique  for  generating  the  phase  screen  does,  in  fact, 

yield  phase  screens  that  have  the  theoretical  phase  structure  function  given 

in  eq.  (8.  6).  The  deviation  at  large  separations  is  caused  by  the  fact  that 

the  total  extent  of  the  phase  screen  used  in  these  studies  is  64  x 1,  106  = 71  cm, 

whereas  the  outer  scale  is  2 m.  Larger  screens  are  needed  to  obtain 

better  agreement  at  large  separation  values.  We  feel  that  this  moderate 

discrepancy  at  large  separations  is  not  important  in  our  simulation  studies, 

because  the  degree  of  correction  obtainable  with  a deformable  or  piston 

plus  tilt  mirror  does  not  depend  significantly  on  the  large  scale  properties 

of  the  turbulence.  (2)  A comparison  on  the  theoretical  curve  for  the  phase 

structure  function  given  in  eq.  (8.  6)  with  the  corresponding  curve  for  the 

infinite  outer  scale  result  reveals  that  the  mean  square  fluctuation  between 

two  points  is  significantly  less  when  L = 2 m than  it  is  when  L = a.  It 

o o 

is  our  contention  that  the  phase  screens  that  we  are  using,  for  which 

L^  = 2 m,  are  adequate  for  the  simulation  studies  because  the  increased 

rms  phase  fluctuation  level  observed  with  L = co  is  due  to  the  additional 

o 

large  scale  variations  that  occur  under  this  condition. 

Each  phase  screen  was  characterized  by  the  initialization  index  of 
the  random  number  generator.  Thus,  phase  screen  3 corresponds  to  the 
phase  screen  generated  by  the  sequence  of  random  numbers  obtained  with  the 
initialization  index  of  3.  Similar  correspondence  exists  for  phase  screen 
32.  The  level  of  turbulence  considered  in  the  simulations  corresponds  to 
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phase  coherence  lengths  of  10  cm  and  5 cm.  An  option  was  incorporated 
into  the  simulation  scheme  which  allowed  more  than  one  COAT  corrected 
image  to  be  added.  In  this  way,  phase  corrected  images  for  many  phase 
screens  could  be  averaged.  Both  single  run  (i.  e. , one  phase  screen)  and 
five  runs  averaged  (i.  e. , averaged  over  five  separate  phase  screens)  images 
were  simtilated.  In  addition,  the  simulations  performed  did  not  include 
the  effect  of  noise  in  the  system. 

The  object  used  in  the  simulations  was  derived  from  a digitized  tape 
of  the  intensity  reflection  of  a scaled  model  of  a satellite,  provided  by 
B.  L.  McGlammery.  The  object  possessed  a dynamic  range  of  1296 
intensity  levels.  Originally,  the  object  was  recorded  in  a 2 56  x 256  mesh- 
point  format.  Because  of  limitations  imposed  by  core  space  in  the  computer, 
the  object  was  converted  to  a 128  x 128  format  by  utilizing  only  that  portion 
of  the  object  which  was  pertinent  to  our  simulation  studies,  A photograph 
of  the  object  used  is  shown  in  Fig.  VIII-2.  The  display  is  limited  to  only 
16  levels  of  grey,  thereby  severely  restricting  the  faithful  reproduction  of 
the  object  because  of  its  large  (1296  levels)  dynamic  range.  In  order  to 
overcome  this  dynamic  range  limitation,  enhancement  schemes  were  applied 
to  the  object  before  display.  Similarly,  the  COAT  reconstructed  images 
were  also  enhanced  prior  to  display. 

The  enhancement  technique,  which  proved  to  be  most  useful,  is 
known  as  the  local  area  bright  less  and  gain  control  (LABGC)  technique. 

This  technique,  developed  by  D.  J.  Ketcham,  tabulates  the  intensity 
distribution  of  the  pixels  within  a localized  window  in  the  image,  calc^llates 
the  statistical  properties  therein,  and  applies  a linear  gain  to  pixel  values 
below  a threshold  value,  determined  by  the  statistics  within  that  window. 
Pixel  values  above  the  threshold  value  are  given  a gain  of  one.  Thus,  the 
low  intensity  features  of  the  image  are  highlighted  while  the  high  intensity 
features  are  unaltered  to  prevent  saturation.  An  illustration  of  the  gain 
curve  is  shown  in  Fig,  VlII-3.  In  the  enhancement  process  the  replacement 
of  the  old  pixel  values  by  the  "enhanced"  values  is  carried  out  for  only  the 
center  portion  of  the  localized  window  in  order  to  avoid  the  mosaic  effect 
which  would  arise  if  the  entire  window  area  is  replaced  by  the  new 
enhanced  values.  The  window  is  then  displaced  by  one  mesh  point,  and 


t 
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the  entire  process  started  over  again  for  the  new  window  until  the  entire 
image  has  thus  been  enhanced  in  this  fashion.  In  our  enhancement  schemes, 
a 9 X 9 window  was  used  ard  a central  2x2  area  was  used  for  the  updating, 

A detailed  discussion  of  the  LABGC  technique  is  behond  the  scope  of  this 
report  and  can  be  found  in  Ref.  11.  The  maximum  gain,  G , , and  the  intcn- 
sity  values  which  determine  in  Fig.  VIII-3  are  input  parameters,  and 
can  be  varied  to  achieve  the  best  enhancement.  It  should  be  emphasized 
that  the  simulation  scheme  has  been  carried  out  for  a noise-free  system. 

In  the  presence  of  noise,  the  enhancement  scheme  may  bring  out  nois%; 
related  features  which  are  not  part  of  the  original  object. 

The  results  of  the  simulation  studies  are  shown  in  Fig.  VIII-4,  with 
the  parameters  of  the  corresponding  simulations  listed  in  Table  VIII- 1. 

Once  chosen,  the  gain  and  threshold  values  used  for  these  enhanced  pictures 
were  kept  constant  for  all  cases.  They  were  chosen,  after  many  trial-euid- 
error  processes,  as  presenting  the  most  pleasing  visual  images  of  the 
GOAT  corr  ected  data,  i.  e.  , the  details  of  the  low  intensity  features  of  the 
image  have  been  properly  highlighted.  It  is  seen  from  Table  Vlll-i  that 
Fig.  VTIi-4  contains  the  images  reconstructed  from  consideration  of  phase 
screen  3.  Figures  VIII-4{a)  and  VIII-4(b)  are  simulations  in  which  the 
phase  coherence  length  is  10  cm;  while  Figs.  VIII-4(c)  and  VIII-4(d)  are  for 
5 cm  phase  coherence  length.  Figures  VIII-4(a)  and  VIII-4(c)  employed 
506  actuators  in  the  simulations,  thereby  representing  a finely  segmented 
mirror;'  Figs.  VIII-4(b)  and  VUI-4(d)  employed  225  actuators  and  represent 
a coarsely  segmented  mirror.  All  simulations  in  Fig.  VIII-4  were  carried 
out  for  only  one  phase  screen.  The  residual  rms  phase  error  for  each  case 
is  shown  in  the  last  column,  and  range  from  0.  582  radians  to  1.  25  radians. 

It  is  seen  that  there  is  a general  improvement  in  image  quality  when  the 
residual  phase  error  is  decreased  significantly.  This  is  shovm  by  com- 
paring Figs,  VIII-4(c)  and  VUI-4(d)  with  Fig.  VIII-4(a).  However,  when 
the  residual  phase  errors  are  not  significantly  different,  the  improvement 
in  image  quality  becomes  less  apparent,  if  at  all.  This  can  be  seen  by 
comparing  Fig.  VIII-4(c)  with  Fig.  VIII-4(d),  or  VIII-4(a)  with  VIII-4fh). 
Similar  remarks  hold  for  the  other  figures  of  Table  I. 
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Table  VIII- 1.  Simvilation  Conditions 


Figure 

Phase 

Screen 

Phase  Coherence 
Length,  cm 

Number  of 
Actuators 

Number  of 
Runs 

Residual 
4*  rms 

VIII-4a 

3 

10 

506 

1 

0.  582 

vm-4b 

3 

10 

225 

1 

0.  750 

VIII-4C 

3 

5 

506 

1 

0.  969 

VIII-4d 

3 

5 

225 

1 

1.  25 

VIII- 5a 

32 

10 

506 

1 

0.  577 

VIII- 5b 

32 

10 

225 

1 

0.790 

vm-5c 

32 

5 

506 

1 

0.  962 

VIII- 5d 

32 

5 

225 

1 

1.  32 

VIII-6a 

3 

10 

506 

5 

0.  590 

VIII-6b 

3 

10 

225 

5 

0.  793 

vni-6c 

3 

5 

506 

5 

0.  984 

VIII-6d 

3 

5 

225 

5 

1.  32 

Vin-7a 

32 

10 

506 

5 

0.607 

VIII-7b 

32 

10 

225 

5 

0.836 

VIII- 7 c 

32 

5 

506 

5 

1.  01 

VIII- 7d 

32 

5 

225 

5 

1.  39 

An  examination  of  the  enhanced  images  of  Figs.  VIII-4  through  7 
shows  that  overall,  when  the  residual  phase  errors  are  significantly 
different  in  magnitude,  there  is  a correlation  between  image  quality  and 
residual  rms  phase  error.  For  example,  the  image  quality  is  significantly 
different  between  the  -a  and  -c  or  -d  images  of  each  figure.  Here,  the 
residual  rms  phase  errors  differed  by  nearly  a factor  of  2 to  2.  5.  How- 
ever, as  the  difference  in  the  phase  errors  is  decreases,  it  becomes  more 
difficult  to  detect  any  obvious  improvements  in  image  quality.  This  point 
is  readily  illustrated  by  Figs.  VIII-4(c)  and  (d),  where  (d)  has  a larger 
value  of  residual  phase  error,  but  it  is  not  obvious  that  the  image  quality 


(c)  (d) 


Fig.  VIII-6.  Degraded  image,  smoothed  using  five  separate 
representations  of  random  phase  screen  No.  3 
for  conditions  in  Table  VIII- 1. 


is  worse  than  Fig.  VIII-4(c).  Other  figures  show  a gradual  degradation  in 

image  quality  as  residual  phase  error  is  increased  slightly.  This  effect 

seems  to  be  more  pronounced  for  lower  values  of  A . The  multiple  images 

' rms  ■ 

which  sometimes  appear  in  these  figures  are  due  to  the  side-lobes  associated 
with  the  segmented  mirror.  As  the  number  of  actuators  is  decreased, 
thereby  necessitating  a larger  block  of  mesh  points  for  each  segment,  the 
period  of  the  segmented  blocks  is  increased.  Thus,  the  sidelobes  appear 
closer  together  and  introduce  closer  spaced  multiple  ’mages  as  in  the 
-b  and  -d  parts  of  each  figures. 
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IX. 


REPORT  CONCLUSIONS 

This  final  report  for  the  predetection  compensated  imaging  study 
provides  the  opportunity  to  review  the  major  conclusions  and  the  signifi- 
cant accompl'shments  of  the  study.  The  nature  of  the  study  itself  changed 
very  early  with  the  contract  redirection  to  emphasize  theoretical  evaluation 
of  the  candidate  predetection  compensated  imaging  concepts. 

One  of  the  most  significant  conclusions  of  the  study  v^as  the 
theoretical  evaluation  of  the  sensitivity  of  the  following  wavefront  error 
sensors: 

• Hartman  test  system 

• Pinhole  interferometer 

• Integral  of  intensity  through  a m ask  (pinhole  image 
quality  sensoi). 

The  Hartman  test  system  was  shown  to  be  a very  sensitive  wavefront 
error  sensor.  Although  its  sensitivity  is  good,  the  pinhole  interferometer 
was  shown  to  have  difficulties  with  wavefront  error  connection  when  viewing 
a multiple  glint  target.  The  pinhole  image  quality  sensor  with  one  detector 
was  found  to  have  extremely  poor  sensitivity. 

However,  another  major  development  in  the  study  was  the  evaluation 
of  the  multiple  aperture  interlinked  (MAIL)  predetection  compensated  imag- 
ing system.  This  MAIL  concept  obtains  a significantly  improved  sensitivity 
compared  with  the  pinhole  image  quality  sensor  by  employing  multiple 
detectors  to  sense  wavefront  error  within  blocks  of  the  full  telescope 
aperture.  The  MAIL  system  was  the  candidate  concept  for  evaluation  in  the 
first  and  second  Standard  Noise  Calculations. 
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APPENDIX 


CALCULATION  I , THE  DEGREE  OF  COHERENCE 
m 


This  appendix  calculates  the  intensity  and  depth  of  modu- 
lation of  one  output  of  an  interferometer  whose  inputs  are  I^ 

and  I . 
w 

The  degree  of  coherence  between  I^  and  I^  can  be  expressed 
by  writing  | 

ej^(t)  = e^(t)  + ej^(t)  I 

and  1 

f 

- ae^it'i  t 62  (t)  i 

where  the  e(t)  are  instantaneous  electric  fields  and  s^,e^, 
and  62  have  no  mutual  coherence: 

<e^  (t)e^  (t) ''  ~ 0 i 

i 

I 

<e^(t)e2(t)>  ~ 0 I 


<eQ{t)e2(t)>  - 0 


The  angular  brackets  denote  time  average  over  about  a 

microsecond  (i.e.,  short  with  respect  to  the  I-Coat  system 

-15 

response,  but  long  compared  v;xth  the  10  ' sec  period  of  light 

oscillations) . The  electric  fields  are  normalized  sc  that 


I^  = <e^(t)  > for  i = n,  w,  o,  1,  2. 


then 


and 


^n  = ^o  ^1 


I = a I^  + I^. 
w o 2 


H?ECESI]^  TPAGS 


The  output  of  an  interferometer  that  combines  I and  I is 

n w 

= <[e^(t)  + e^(t+A)]^> 

= ^n  ^ ^ 2<e^(t)e^(t+A)> 

= I + I + 2a  e (t)e  (t+A)>  + equal  to 
n w o o 

zero,  where  A is  the  path  difference  divided  by  the  speed  of 
light.  Assuming  an  average  light  frequency  v,  and  some  degree 
of  chromatism,  for  small  A we  have 

<e^ (t) e^ (t+A) > ~ cos  Av, 

then 

= 1/2  (I^  ^“^o  * 

The  average  intensity  is  + I^,  as  expected,  and  2al^  is 

the  peak  modulation.  We  define  = 2al 

mo. 

We  now  calculate  and  evaluate  I for  various  cases.  A 

o 

cross  section  of  the  intensity  distribution  before  and  after 
the  narrow  and  wide  pinholes  is  shown  to  the  right  for  each 
case.  Since  both  the  narrow  and  wide  pinholes  are  approxi- 
mately unresolvable  by  each  subaperture,  each  of  the  n inter- 
ferometers has  as  its  two  inputs  1/n  of  the  total  light  pene- 
trating the  wide  and  narrow  pinholes. 

1.  Cases  I and  II,  Small  c 

Since  all  the  light  originates  from  an  unresolved 

point  (small  o),  I and  I are  fully  coherent: 

n w 


176 


0 


4845-13 


^2 


I = I 
o n 


a = .^  /I'n 
w n 


so 


1=2  /nr 

m w n 


2 . Cases  I and  II » Large  o 

Since  the  system  is  converged,  different  areas  of  the 
image  on  the  pinholes  originate  from  different  areas  on  the 
target  and  are  thus  uncorrelated.  So  the  shaded  portion  of 
che  image  passing  the  wide  pinhole  is  coherent  with  the  light 
passing  the  narrow  pinhole: 


I = I 
o n 


a2. 


V 


2 


so 


2 ^ 


I = ^ /I  I 
m V w n 


4845-IZ 


3 . Cases  III  and  IV,  Small  g 

As  in  1.  above,  all  the  light  originates  from  an 

unresolved  point  so  I and  I are  fully  coherent: 

n w 


177 


4845-11 


1=2  /Tir 

m w n 


_/V^  A 


4 . Cases  III  and  IV,  Large  o 

Here,  the  unconverged  condition  causes  the  light  from 

* 

each  target  point  to  be  spread  to  a width  Vd^.  The  light 
penetrating  the  narrow  pinhole  is  therefore  a (roughly)  uni- 
form mixture  of  light  from  a circular  region  on  the  target 
of  angular  diameter  Vd^^/H.  The  light  penetrating  the  wide 
pinhole  is  a mixture  of  light  from  a region  on  the  target  of 
angular  diameter  2Vd^/£,  with  the  mixture  weighted  by  a radial 
function  which  gives  the  overlap  of  two  equal  circles  of 
angular  diameter  Vd^/2,.  One  of  these  circles  is  the  spread 
caused  by  nonconvergence  and  the  other  is  the  wide  pinhole. 

The  weighting  functions  for  the  narrow  and  wide  pinholes  look 
like 


4845-14 


NARROW 

PINHOLE 


WIDE 

PINHOLE 


-Vd„/I 


RADIUS,  rodions  ■ 
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I =:  <^6  > 

W W 


,(r)  dr 


The  interferometer  output  is  then 


I = I + I + T cos  (a-U)  , 

c n w m V / f 


I H 2 /f  (r)  f (r)  I . 

m / n-  r r 


Recalling  our  original  definition  of  the  target,  we  have 


a e 


I (2  /o  2 

- r /2o 


If  0 is  large,  then  is  nearly  constant  when  f^^  and  f^  are 
nonzero,  and 


1=21 

m 


00 

r=0  [ ^n^^ 


) f„(r)  d'-r  , 
w - 


I = I 
n r 


00 

=o/ 


) dr  , 


00 

- -"r^O  / 


f (r)  d^r  . 
w - 


Therefore , 


1=2  /I  I 
m n w 


00 

/ v/f  (r)  f I 
j n ~ w 


~-r  ^2 
r)  d r 


If 


For  our  functions  f and  f 

n w 


I - 1.4  /I i 
m n w . 
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